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C3 ; ABSTRACT 

. We present the final results from our deep HST imaging study of the host galaxies of 

' radio-quiet quasars (RQQs), radio- loud quasars (RLQs) and radio galaxies (RGs). We 

! describe and analyze new WFPC2 i?-band observations for 14 objects which, when 

combined with the first tranche of HST imaging reported in McLure et al (1999), 
(N ■ provide a complete and consistent set of deep, red, line-free images for statistically- 

matched samples of 13 RQQs, 10 RLQs and 10 RGs in the redshift band 0.1 < z < 0.25. 
We also report the results of new deep VLA imaging which has yielded a 5 GHz 
detection of all but one of the 33 AGN in our sample. 
. Careful modelling of our images, aided by a high dynamic-range point-spread 

' function, has allowed us to determine accurately the morphology, luminosity, scale- 

^-H i length and axial ratio of every host galaxy in our sample. Armed with this information 

■ we have undertaken a detailed comparison of the properties of the hosts of these 3 
' types of powerful AGN, both internally, and with the galaxy population in general. 

We find that spheroidal hosts become more prevalent with increasing nuclear lu- 
I ' minosity such that, for nuclear luminosities My < —23.5, the hosts of both radio-loud 

O ' and radio-quiet AGN are virtually all massive ellipticals. Moreover, we demonstrate 

^ , that the basic properties of these hosts are indistinguishable from those of quiescent, 

^ ■ evolved, low-redshift ellipticals of comparable mass. This result rules out the possibil- 

I ity that radio-loudness is determined by host-galaxy morphology, and also sets severe 

^ . constraints on evolutionary schemes which attempt to link low-z ULIRGs with RQQs. 

' Instead, we show that our results are as expected given the relationship between 

I black-hole and spheroid mass established for nearby galaxies, and apply this relation 

■ to estimate the mass of the black hole in each object. The results agree remarkably 
■ " " ' well with completely-independent estimates based on nuclear emission-line widths; all 

the quasars in our sample have Mi,h > 5 x 10* Mq while the radio-loud objects are 
confined to Mth > 10^ M©. This apparent mass-threshold difference, which provides 
a natural explanation for why RQQs outnumber RLQs by a factor of 10, appears 
to refiect the existence of a minimum and maximum level of black-hole radio output 
which is a strong function of black- hole mass (oc M^f^'^'^). Finally, we use our results 
to estimate the fraction of massive spheroids/black- holes which produce quasar-level 
activity. This fraction is ~ 0.1% at the present day, rising to > 10% at z ~ 2 — 3. 

Key words: galaxies: active - galaxies: photometry - infrared: galaxies - quasars: 
general - black hole physics 



1 INTRODUCTION 

Studies of the host galaxies of low-redshift quasars can en- 
able us to define the subset of the present-day galaxy popu- 
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lation which is capable of producing quasar-level nuclear ac- 
tivity. This is of obvious importance for constraining phys- 
ical models of quasar evolution (Small & Blandford 1992; 
Haehnelt & Rees 1993; Kauffman & Haehnelt 2000), and 
for exploring the connection between black-hole and galaxy 
formation (Silk & Rees 1998, Fabian 1999, Franceschini et 
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al. 1999, Granato ot al. 2001, Kormendy & Gcbhardt 2001). 
Such observations arc also of value for testing unified models 
of radio-loud AGN (e.g. Peacock 1987, Barthel 1989, Urry 
& Padovani 1995), constraining possible evolutionary links 
between ULIRGs and quasars (Sanders & Mirabel 1996), 
exploring the origin of radio-loudness (Blandford 2000), and 
as a means to estimate the masses of the central black holes 
which power the active nuclei (McLure et al. 1999). 

Our view of low-redshift quasar hosts has been clari- 
fied enormously over the last five years, primarily due to 
the angular resolution and dynamic range offered by the 
Hubble Space Telescope (HST). After some initial confu- 
sion, recent HST-based studies have now reached agreement 
that the hosts of all luminous quasars {My < —23.5) are 
bright galaxies with L > L* (Bahcall et al. 1997, McLure 
et al. 1999, McLeod & McLeod 2001). However, it can be 
argued, (with considerable justification) that this much had 
already been established from earlier ground-based studies 
(e.g. Smith et al. 1986, Veron-Cetty & Woltjer 1990, Taylor 
et al. 1996). 

In fact, as first convincingly demonstrated by Disney et 
al. (1995), the major advance offered by the HST for the 
study of quasar hosts is that it allows host galaxies to bo 
mapped out over sufficient angular and dynamic range for a 
de Vaucouleurs r^''*-law spheroidal component to be clearly 
distinguished from an exponential disc, at least for redshifts 
z < 0.5. This is not to suggest that AGN host-galaxy mor- 
phological discrimination has proved impossible from the 
ground. Indeed for lower- luminosity AGN aX z < 0.1, such as 
Seyfert galaxies, ground-based imaging has proved perfectly 
adequate for this task (e.g. Hunt et al. 1999) and in fact 
some early ground-based attempts to determine the mor- 
phology of low-redshift quasar hosts have also proved to be 
robust (e.g. Smith et al. 1986). However, to ensure an unbi- 
assed comparison of RQQ and RLQ hosts it is necessary to 
study host galaxies at « > 0.15 and to be able to determine 
host-galaxy morphologies for quasars with luminosities up 
to Mv < —26. Even by moving to the infrared to minimize 
nuclearihost ratio, Taylor et al. (1996) found that this could 
not be reliably achieved with typical ground-based seeing. 

Nevertheless, groat care needs to be taken to extract the 
full benefit of HST imaging of quasar hosts. In particular, 
deep observations are required to detect the extended low 
surface-brightness emission of even a massive host galaxy at 
z ~ 0.2 to a radius of several arcsec from the nucleus. Unfor- 
tunately however, this inevitably leads to saturation of the 
nucleus, making accurate characterization of the luminos- 
ity of the central source impossible. This is crucial because, 
at the depths of interest for reliable host-galaxy character- 
ization, scattered light in the WFPC2 PSF still makes a 
significant contribution to surface brightness out to an an- 
gular radius ~ 10 arcsec (McLure, Dunlop & Kukula 2000). 
As demonstrated by McLeod & Rieke (1995), these prob- 
lems of surface brightness bias, saturation, and inadequate 
knowledge of the large-angle properties of the true WFPC2 
PSF, can explain much of the confusion produced by the 
first studies of quasar hosts undertaken after the correction 
of the HST optics with COSTAR (e.g. Bahcall, Kirhakos & 
Schneider 1994). 

In this paper we present the final results from our 34- 
orbit Cycle-6 imaging study of quasar hosts, which was care- 
fully designed to avoid these problems. Specifically, we ac- 



quired images of each quasar spanning a wide range of inte- 
gration times (to allow an unsaturated, high dynamic-range 
image of each object to be constructed) and devoted an en- 
tire orbit to the construction of the necessary high dynamic- 
range PSF (via observations of a star of similar colour to the 
quasar nuclei, imaged at the same location on the same WF 
chip). Results from the first half of this programme were 
reported in McLure et al. (1999), where images for 19 ob- 
jects from our 33-source sample were presented, modelled 
and analyzed. Here wc present and model the images for 
the 14 targets which were observed in the latter half of 1998 
and in 1999, and then summarize and discuss the results de- 
rived from the analysis of the completed sample. The results 
presented in this paper thus complete, extend and in several 
cases supercede those presented in McLure et al. (1999) (e.g. 
estimated black-hole masses for all objects are now calcu- 
lated using more recent estimates of the black-hole:spheroid 
mass relation, yielding significantly lower values than were 
calculated by McLure et al. based on the relation presented 
by Magorrian et al. (1998)). 

Several other substantial studies of low-redshift quasar 
hosts have now been undertaken with the HST (e.g. Bahcall, 
Kirkhados & Schneider 1997; Hooper, Impey & Foltz 1997; 
Boyce et al. 1998, McLeod & McLeod 2001). However, one 
unique feature of the present study is the deliberate focus 
on a comparison of the hosts of the three main classes of 
powerful AGN, namely radio-quiet quasars (RQQs), radio- 
loud quasars (RLQs) and radio galaxies (RGs). Moreover, 
we have ensured that this comparison can be performed in 
an unbiassed manner by confining our sample to a narrow 
range in redshift (0.1 < z < 0.25) and requiring that the in- 
dividual sub-samples are matched in terms of their luminos- 
ity distributions (optical luminosity in the case of the RQQ 
and RLQ sub-samples, and radio luminosity in the case of 
the RLQ and RG sub-samples - see Dunlop et al. (1993), 
McLure et al. (1999) and Section 2 for further details). An- 
other strength of this study is the wealth of pre-existing 
data at other wavelengths, as detailed in Section 2. This 
has allowed us to maximise the impact of the HST imaging 
(through, for example, the determination of R — K colours 
for all the host galaxies). Finally it is worth emphasizing 
that in this study we have sought to extract the properties of 
the stellar population which dominates the mass of the host 
galaxy. Thus, while we do include a statistical analysis of the 
prevalence of peculiar features such as tidal tails, we have 
endevoured to minimize the distorting effect of the more 
transient activity by insisting on line-free imaging longward 
of the 4000A break, and masking out obvious assymmetries 
prior to modelling the host morphology. This approach, cou- 
pled with the careful design of our observations, is the most 
likely explanation for why the results presented in this pa- 
per are generally cleaner, and more homogenous than the 
results of many other recent studies. 

The layout of this paper is as follows. In Section 2 we re- 
view the main properties of the matched RG, RLQ and RQQ 
samples, and summarize the wealth of supporting ground- 
based multi-frequency data which now exists for these ob- 
jects. In Section 3 we give details of the HST observations 
and briefly review the process of data reduction and PSF de- 
termination (McLure et al. 1999; McLure, Dunlop & Kukula 
2000). Then, in Section 4 we give a brief description of the 
new VLA observations of the RQQs in our sample which es- 
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caped previous radio detection by Kukula et al. (1998), and 
present new 5 GHz flux densities and positions. In Section 5 
we present the new HST data and briefly summarize the ap- 
proach taken to modeUing the WFPC2 images; the images 
and 2-dimonsionaI model fits for the 14 now objects are pro- 
vided in Appendix A (along with brief notes on each source), 
with observed and fitted luminosity profiles presented in Ap- 
pendix B. Next, in Section 6 we analyze the results of mod- 
elling the HST images of the complete sample, and assess 
the implications of our results for unified models of AGN. 
In Section 7 we explore how quasar hosts relate to the gen- 
eral population of massive galaxies, and demonstrate via the 
Kormcndy relation that the quasar hosts arc drawn from 
the population of massive 'boxy' ellipticals which arc found 
predominantly in cluster environments. We also include a 
comparison of the environments and interaction statistics 
for quasar hosts and brightest cluster galaxies. In Section 
8 we exploit our results to investigate the properties of the 
central engines which power these luminous AGN, deriving 
estimates of black-hole mass and Eddington ratio, and ex- 
ploring possible clues to the origin of radio-loudness. Finally, 
in Section 9 we briefly discuss the implications of our results 
for the cosmological evolution of quasars and massive galax- 
ies. Our conclusions are summarized in Section 10. 

For ease of comparison with previous work we adopt an 
Einstein-de Sitter universe with Ho = 50 kms~^Mpc~^ for 
the calculation of physical quantities throughout this paper. 



2 SAMPLE AND ASSOCIATED 
OBSERVATIONS 

The HST imaging observations reported here complete the 
imaging of the full sample of 33 objects (10 RLQs, 13 RQQs 
and 10 RGs) defined for this study. This sample was se- 
lected from the slightly larger (40-source) sample imaged in 
the infrared by Dunlop et al. (1993) and Taylor et al. (1996) 
through the imposition of the slightly more restrictive red- 
shift limits 0.1 < 2; < 0.25. As described in McLure et al. 
(1999), this restriction in redshift range ensures that our _R- 
band imaging through the F675W filter is not contaminated 
by the presence of strong emission lines such as [Olll], or 
Ha. The main effect of this additional redshift restriction 
is to exclude a small number of objects from the original 
sample of Dunlop et al. (1993) which have 0.25 < z < 0.35. 
However, this has not significantly compromised the origi- 
nal statistical merits of this sample, namely that the RLQ 
and RQQ sub-samples are matched in terms of optical lumi- 
nosity, and that the RLQ and RG samples are matched in 
terms of radio luminosity and radio spectral index (Dunlop 
et al. 1993). 

Of the 33 sources in this sample, 19 were observed dur- 
ing the first year of our Cycle 6 allocation. The observations 
and analysis of these objects were presented by McLure et 
al. (1999). The remaining 14 objects for which the observa- 
tions are presented and analyzed in this paper are listed in 
Table 1, along with the dates on which they were observed 
with the HST. 

Also included in this paper are new, deep VLA obser- 
vations of a subset of our RQQ sample. These 4.8 GHz ob- 
servations, the results of which are presented in Section 4, 
go a factor of 3 deeper than the observations of Kukula et 



Table 1. Observing dates for the objects presented in this paper. 
Note that, despite its archive designation as 3C59, 0204-1-292 is 
in fact now classified as an RQQ (see Taylor et al. 1996, and the 
radio luminosity quoted in Table 3). 



Object 


HST Archive 


Type 


Observing 




designation 




date 


0307+169 


3C79 


RG 


Jul 10 1998 


0230-027 


PKS0230-027 


RG 


Sep 25 1998 


1342-016 


1342-016 


RG 


Dec 10 1998 


1215-1-013 


1215+013 


RG 


Jan 02 1999 


1215-033 


1215-033 


RG 


Jan 06 1999 


1330+022 


1330+022 


RG 


Apr 13 1999 


1217+023 


PKS1217+02 


RLQ 


Jul 07 1998 


1020-103 


PKS1020-103 


RLQ 


Jul 12 1998 


2135-147 


PKS2135-14 


RLQ 


Oct 19 1998 


2355-082 


PKS2355-082 


RLQ 


Oct 19 1998 


0204+292 


3059 


RQQ 


Jul 09 1998 


1549+203 


1E15498+203 


RQQ 


Sep 03 1998 


2215-037 


EX2215-037 


RQQ 


Sep 26 1998 


0052+251 


PG0052+251 


RQQ 


Nov 06 1998 



al. (1998) which were utilised by McLure et al. (1999). The 
important outcome of these observations is that we have 
now detected all but one of the RQQs in this sample at ra- 
dio wavelengths. Those new radio detections, coupled with 
completion of the HST observations, have allowed us to re- 
investigate and clarify a number of the relations between 
optical and radio properties which could only be tentatively 
explored by McLure et al. (1999). 

We have also obtained improved infrared (UKIRT K- 
band) images for a small number of the more luminous 
quasars in our sample since the publication of McLure et 
al. (1999). These have boon published in McLure, Dunlop 
& Kukula (2000), but the results of modelling these images 
arc utilised in this paper to assist in the improved analysis 
of galaxy colours which is presented in Section 6.6. 

Finally we note that H/3 emission-line spectroscopy of 
the RLQ and RQQ samples discussed here has recently been 
completed by McLure & Dunlop (2001). The results of this 
spectroscopic study are referred to in the discussion of black- 
hole mass estimation presented in Section 8. 



3 HST OBSERVATIONS 

The observations were made with the Wide Field & Plane- 
tary Camera 2 ( WFPC2; Trauger et al. 1994) on the Hubble 
Space Telescope through the F675W filter. The filter spans a 
wavelength range of 877A from 6275.5 to 7152. 5A roughly 
equivalent to standard R band. This filter was selected in 
preference to a wider filter because it allowed both [Oiii] 
and Ha emission lines to be excluded from the band-pass 
for source redshifts in the range 0.1 < 2: < 0.25, and thus 
ensured that a clean measure could be made of the level 
of continuum light emitted by the quasar host galaxy at 
wavelengths longwards of the 4000A break. As in McLure 
et al. (1999) target sources were centred on the WF2 chip, 
which was chosen in preference to WFl or WF3 because of 
its marginally superior performance over the period imme- 
diatly prior to our observations. 

Deep sensitive images of the host galax:ies are obviously 
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desirable, but the necessary long, background-limited ex- 
posures inevitably mean that, in the case of the quasars, 
the central source becomes saturated. Prom such saturated 
images it is extremely difficult to reliably disentangle host- 
galaxy omission from the contribution of the PSF-convolved 
nuclear source. 

Slightly different strategies were therefore adopted for 
the quasar and radio-galaxy samples. For the quasars, expo- 
sures of 5, 26 and 3x600 seconds were taken. The short ex- 
posures ensured that at least one unsaturated image of each 
quasar would be obtained, thus enabling an accurate mea- 
surement of the central flux density. The three 600-second 
exposures each yielded a 3a surface-brightness sensitivity 
of fiR = 23.8 mag. arcsec^'^ px—1, and their comparison fa- 
cilitated reliable cosmic ray removal using standard IRAF 
tasks. With azimuthal averaging, the combined 1800-second 
image of each quasar allows extended emission to be traced 
reliably down to a surface- brightness level of /lih > 26 mag. 
arcsec"^. 

For the radio galaxies there was little danger of sat- 
uration and so short exposures were not required. Three 
700-second exposures were therefore obtained for each ra- 
dio galaxy. Any remaining time in the orbit was filled with 
a shorter exposure of flexible length (usually 40 to 100 
seconds). Calibration was carried out using the standard 
pipeline. 

As described in McLure et al. (1999) we devoted one 
orbit of our allotted HST time to constructing a deep, un- 
saturated stellar PSF using the F657W filter, with the star 
centred on exactly the same part of the WF2 chip as the tar- 
get objects. Full details (magnitude, colour etc) of the chosen 
star can be found in McLure et al. (1999), along with a de- 
tailed description of how we constructed a properly-sampled 
PSF of the required large dynamic range using dithered ob- 
servations of a series of exposures ranging from 0.23 to 160 
seconds. 

Here we simply re-emphasize that the resulting high 
dynamic range PSF, tailored as closely as possible to match 
our quasar observations, has proved to be absolutely crucial 
in allowing the reliable extraction of host-galaxy parameters 
from the WFPC2 images. 



4 NEW RADIO OBSERVATIONS OF RQQS 

Kukula et al. (1998) reported radio continuum observations 
of a sample of 27 low-redshift radio-quiet quasars, including 
all of the RQQs in the current sample, made with the Very 
Large Array (VLA) at 1.4, 4.8 and 8.4 GHz. Although flux 
densities were obtained for the majority of the sample, six 
objects remained undetected at all throe frequencies (the 
flux-density limit of the survey was ~ 0.2 mjy at 4.8 GHz). 

To remedy this situation additional high-sensitivity ob- 
servations were made of these six RQQs in August 1999 
(05/08/99) at 4.8 GHz, again using the VLA in its high- 
resolution 'A' configuration. Data reduction was carried out 
using the standard procedure within AIPS. The angular res- 
olution of the maps is ~ 0.4" (FWHM) at this frequency, 
corresponding to a physical size of ^ 1.5 kpc at the rcdshift 
of the quasars - thus neatly encompassing the nuclear region 
of the quasar host. 

The results of these observations are listed in Table 2. 



Radio sources are detected within ~ 1" of all bar one of 
the RQQs (positional accuracies are estimated to be within 
~ 100 mas). Due to scheduling problems 1549-1-203 could 
only be allocated half the time given to the other five RQQs 
and the map of this object suffers from a correspondingly 
higher noise level. 



5 NEW HST RESULTS 

The WFPC2 F675W images, two-dimensional model fits, 
and the model-substracted residual images of the 14 new 
objects are presented in Appendix A in Figs Al A14, along 
with brief notes on each individual source. Comparable im- 
ages and notes for the other 19 objects in the sample can 
be found in McLure et al. (1999). The observed luminosity 
profiles for the 14 new objects are presented in Appendix B 
in Figs B1-B14, along with the best-fitting model profiles 
extracted from the two-dimensional model fits. 

Full details of the two-dimensional model procedure 
which we have used to determine the properties of the host 
galaxies can be found in McLure, Dunlop & Kukula (2000) , 
along with the results of extensive tests of its ability to re- 
claim host-galaxy parameters from simulated data based on 
a wide range of host-galaxy:nucleus combinations at differ- 
ent redshifts. As emphasized in McLure, Dunlop & Kukula 
(2000), the success of this modelling depends on an accu- 
rate high dynamic range PSF, and the construction of an 
accurate error frame for each quasar image. 

In brief, the modelling of the HST images was carried 
out in three separate stages. The first stage involves assess- 
ing how well the data can be reproduced assuming that the 
host galaxy is either an elliptical galaxy (with a surface- 
brightness distribution described by a de Vaucouleurs r^^*- 
law) or a pure exponential disc. The remaining five param- 
eters (host-galaxy position angle, host-galaxy axial ratio, 
host-galaxy scale-length, host-galaxy luminosity, and nu- 
clear luminosity) are then adjusted until, when convolved 
with the PSF, the model best fits the data as determined 
by x-squared minimization (note that it is not assumed a 
priori that the radio galaxies have a negligible nuclear com- 
ponent). Then, if one assumed galaxy morphology yields a 
significantly better fit than the other, we can say that the 
galaxy is better described by a de Vaucouleurs law or by an 
exponential disc. As with all the modelling performed on 
the HST sample, once the minimum x solution had been 
found, the modelling code weis repeatedly re-started from 
close to the minimum solution, in order to ensure that 
the solution was stable. 

The results of applying this procedure to the new HST 
images are given in Table 3, alongside the results already 
determined by McLure et al. (1999). The striking feature 
of these results, now confirmed with the complete sam- 
ple, is that with the exception of 3 RQQs (0052-f251 and 
the two lowest-luminosity RQQs 0257-1-024 and 2344-1-184) 
every single host galaxy is better described by a de Vau- 
couleurs law. 

In our second approach we have removed the require- 
ment of assuming that the host galaxy can be described as 
either a pure r^^''-law or exponential disc, and allow a sixth 
parameter /3 (where the luminosity profile of the galaxy is 
given by /(r) oc exp{—r^)) to vary continuously. Thus /3 = 1 
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Table 2. High-sensitivity VLA observations of the six previously-undetected RQQs in our sample. Optical positions were measured from 
Digitised Sky Survey plates using the Stsdas package in irap. All radio observations were made in C-band (4.8 GHz) with the VLA in 
A-configuration (angular resolution ~ 0.4") on August 5 1999. The uncertainties in the measured flux densities are given as 3 times the 
rms noise in the image. Radio positions are accurate to within 100 mas. 



Optical Position (J2000) Radio Position (J2000) 4.8 GHz Flux log(P4.8Gffz/ 



Quasar 


z 


My 


RA {h m s) 


Dec (° 


"') 


RA {h m s) 


Dec (° 


"') 


Density /mjy 


/W Hz-isr 


0244+194 


0.176 


-23.55 


02 


47 


40.85 


+19 


40 


57.8 


02 


47 


40.84 


+19 


40 


57.8 


0.18 


± 


0.06 


21.3 


PG 0923+201 


0.190 


-24.56 


09 


25 


54.71 


+19 


54 


04.4 


09 


25 


54.74 


+19 


54 


05.0 


0.14 


± 


0.06 


21.3 


PG 0953+414 


0.239 


-25.36 


09 


56 


52.35 


+41 


15 


22.5 


09 


56 


52.39 


+41 


15 


22.2 


0.25 


± 


0.07 


21.7 


1549+203 


0.250 


-24.51 


15 


52 


02.36 


+20 


14 


00.5 














< 


0. 


12 


<21.4 


2215-037 


0.241 


-23.73 


22 


17 


47.77 


-03 


32 


38.8 


22 


17 


47.72 


-03 


32 


38.5 


0.13 


± 


0.08 


21.4 


2344+184 


0.138 


-23.76 


23 


47 


25.71 


+18 


44 


50.8 


23 


47 


25.77 


+18 


44 


50.7 


0.19 


± 


0.08 


21.1 



should result if the galaxy is best described by a pure expo- 
nential disc, and B = 0.25 should result if the galaxy really 
does follow a pure de Vaucouleurs law, but all values of /3 
are available to the program if this results in an improved 
quality of fit. The results of applying this procedure to the 
HST images are given in Table 4 (again for the complete 
sample) and illustrated in Fig 1. These results are discussed 
in more detail in Section 6.1. 

Finally, an examination of Table 4 and Fig 1 reveals that 
whereas a very clean preference for /3 = 0.25 is displayed by 
the RGs and RLQ hosts, a few of the RQQ hosts (in par- 
ticular, as mentioned above, 0052+251 and the two lowest- 
luminosity RQQs 0257+024 and 2344+184) have best-fit /3 
values which are intermediate between the values of 0.25 
or 1.0 expected for pure elliptical or disc hosts. For this 
reason it was decided that the RQQ hosts should all be 
re-modelled with a 9 free-parameter fit, which allowed for 
the combination of both disc and bulge contributions to 
the host's surface-brightness distribution. For nine out of 
the 13 RQQs this extra freedom still resulted in no signifi- 
cant disc component. However, for four objects (the above- 
mentioned three RQQs plus 0157+001) this procedure pro- 
duced a significantly-improved model fit, and it is the Lhost 
and Lnuc values from these combined fits which are used in 
all of the subsequent analyses. These four bulge-disc com- 
binations are also noted in Table 3, and in that table (and 
in subsequent analysis) it is the scale-length, axial ratio and 
position angle of the dominant component which is adopted 
for these combined-fit objects. 

The luminosity profiles, extracted from the two- 
dimensional model fits are presented for the 14 new objects 
in Appendix B. The profiles are followed out to a radius of 
10", which is representative of the typical outer radii used 
in the modelling ((r) = 11"). 



6 ANALYSIS OF THE FULL SAMPLE 
6.1 Host- Galaxy Morphologies 

Examination of the results presented in Tables 3 and 4 con- 
firms that the somewhat complex observing strategy out- 
lined in Section 3 has successfully allowed the determination 
of host-galaxy morphology for all 33 objects in the sam- 
ple. Also immediately apparent from Table 3 is that the 
huge amount of information available to the modelling code 



has not only allowed a clear morphological preference to be 
made, but can formally exclude the alternative host in all 
cases (Ax^ = 25.7 corresponds to a 99.99% confidence level 
for a 5 parameter fit). 

The results from the modelling of the RG and RLQ 
sub-samples are in good agreement with orientation-based 
unification, with all 20 objects found to have elliptical host 
galaxies. A perhaps more striking feature of these results 
is the extent to which the classic r^''^ de Vaucouleurs law 
provides a near-perfect description of the host galaxies of the 
radio- loud objects (see Fig 1). As can be seen from Table 4 
the best-fitting fH values for the combined RG and RLQ sub- 
samples all lie in the narrow range 0.19 < l3 < 0.26, with 
the RG sample alone displaying an even narrower spread, 
0.21 < f3 < 0.25. This conclusion is further strengthened 
by a comparison of the beta histogram for the three sub- 
samples with that from the /3-modelling tests performed by 
McLure, Dunlop & Kukula (2000). The application of the 
Kolmogorov-Smirnov test to the (3 distribution of the 29 
objects which were found to have single-component elliptical 
host galaxies, and those resulting from the /3-model testing, 
returns a probability of p = 0.23. This confirms that, as far 
as the modelling code is concerned, any differences between 
these hosts and pure elliptical galaxies are not statistically 
significant. The best match between the test results and the 
actual data is for the RG sub-sample (p = 0.39), perhaps 
as would be expected considering the lack of a dominant 
point-source contribution. 

Considering the long-standing belief that RQQs are of- 
ten located in disc galaxies, the results from modelling of 
the RQQ sub-sample are quite unambiguous, with 9 of the 
13 host galaxies showing no evidence for any disc compo- 
nent. Of the four objects which are best matched by a com- 
bined disc/bulge model, the luminosities of 0157+001 and 
0052+251 are dominated by their bulge components, which 
respectively account for 83% and 71% of the total host lu- 
minosity. This means that the number of bulge-dominated 
RQQ hosts is 11 out of the 13 objects. 

It is interesting to note that the two most disc- 
dominated galaxies, 0257+024 and 2344+184, are the hosts 
of by far the lowest luminosity AGN in the entire 33- 
objcct sample. In fact, converting their total luminosity 
(host+nucleus) from the model fits for these two objects 
to the equivalent absolute F-magnitudes gives My = —22.6 
and Mv = —21.9 respectively (assuming V — R = 0.8 at 
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Table 3. The outcome of attempting to model the AGN host galaxies as either an exponential disc or a de Vaucouleurs spheroid. Source 
name and redshift are given in the first two columns, with the logarithm of radio luminosity given in column 3. The preferred host-galaxy 
morphology is given in column 4, with the Ax^ between the chosen and alternative model listed in column 5. In column 6 ri/2 is given 
irrespective of the chosen host morphology. Column 7 lists units of R mag arcsec"^. Columns 8 and 9 list the integrated apparent 

magnitudes of the host galaxy and fitted nuclear component converted from F675W to Cousins R-band, while column 10 gives the ratio 
of integrated galaxy and nuclear luminosities. Columns 11 and 12 give the axial ratio and position angle (east of north) of the best-fiting 
host-galaxy model. 
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z = 0, Fukugita at al. 1995). Given that for this study 
the adopted quasar /Seyfert borderline is Mv ~ —23.0, it 
is clear that these two objects are not actually bona fide 
RQQs. The clear implication from this result is that all 
true quasars, with My < —23.0, reside in luminous bulge- 
dominated hosts, irrespective of their radio power. 

The morphological determinations for the RQQ host 
galaxies have placed on a firm footing the suggestion made 
previously by Taylor at al. (1996) and McLeod & Rieke 
(1995) that the probability of a RQQ having an early-type 
host was an increasing function of the quasar luminosity. Un- 
like the previous ground-based studios, the high resolution 
and temporally stable PSF offered by HST has permitted 
the confirmation of what were hitherto necessarily tentative 
conclusions, due to the uncertainties introduced by ground- 



based seeing conditions. Therefore, a strong conclusion from 
the morphological analysis of from the new HST images is 
that the radio luminosity of an AGN is not directly related 
to host-galaxy morphology. The relationship between host 
morphology and AGN luminosity is explored further in Sec- 
tion 7. 



6.2 Host-Galaxy and AGN Luminosities 

The host and nuclear luminosities axe presented in the form 
of integrated absolute Cousins J?-band magnitudes in Ta- 
ble 5, and in Figs 2 and 3. These have been derived from 
the apparent magnitudes listed in Table 3 which have been 
calculated by integrating the best-fit model components to 
infinite radius, and adopting the F675W fiight-system zero- 
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Table 4. The outcome of the variable- /3 modelling. Column 2 lists 
the host morphology of the best fitting 'fixed model (results 
of which are given in Table 3). The best- fitting values for the /? 
profile parameter are given in column 3. Column 4 gives Ax^, 
which quantifies the improvement in fit off'ered by the variable-/? 
model over that already achieved with the best-fitting pure disc 
or elliptical model. Simulations indicate that /3 can be reclaimed 
to within a typical uncertainty of 0.01 — 0.02. 
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point given by Holtzman et al. (1995). The similarity be- 
tween the F675W filter and the standard Cousins 7?-band 
is such that, with the use of this zero-point, the difference 
is of the order ±0.05 magnitudes. The conversion from ap- 
parent to absolute magnitudes has been performed using k- 
corrcctions assuming spectral indices of a = 1.5 and a = 0.2 
for the host and nucleus respectively {f^ oc i^^"). 

For all but 4 of the sources, the values of M_R(host) 
Usted in Table 5 can be regarded as equivalent to the value of 
Mj{ (bulge). For the 4 objects which benefitted from a com- 
bined disc+bulgc model the best-fitting values of M_r (bulge) 
are as follows: 0052+251: Mij(bulge) = -22.95, 0157+001: 
Mj{(bulge) = -24.32, 0257+024: Mj{(bulge) = -21.48, 
2344+184: MH(bulge) = -22.44. 
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Figure 1. Sub-sample histograms of the best-fit /3 values from 
the variable-^ modelling. The dotted line lies at /3 = 0.25, corre- 
sponding to a perfect de Vaucouleurs law. 



6.2.1 Host-Galaxy Luminosities 

The mean and median integrated absolute magnitudes of 
the best-fit host galaxies in each sub-sample are: 



(Mr) = 

{Mr) = 

(Mr) = 

(Mr) = 



-23.53 ± 0.09 
-23.66 ±0.16 
-23.73 ±0.10 
-23.28 ±0.15 



median=— 23.52 

median=— 23.63 
median=— 23.67 
median=— 23.30 



(ALL) 

(RG) 
(RLQ) 

(RQQ) 



Two features of these results merit comment. First, the 
agreement between the absolute magnitudes of the host 
galaxies of the RG and RLQ sub-samples can be seen to be 
extremely good, with the median figures differing by only 
0.04 magnitudes. This can be interpreted as strong evidence 
in favour of orientation-based radio-loud unification (see also 
Section 6.7). The second obvious feature of these results is 
that these new HST images appear to confirm the tradi- 
tional finding that the hosts of RQQs arc less luminous than 
those of RLQs, although the median difference of 0.37 mag- 
nitudes is a factor of two smaller than the difference typically 
claimed (e.g. Kirhakos et al. 1999). This does not seem to be 
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Table 5. Absolute magnitudes (Mn), and optical-infrared {R — K) colours of the best-fitting host galaxy and nuclear component for 
each AGN. Column 2 gives the K-band absolute magnitudes (M^) of the total host galaxy derived from the current modelling of the 
HST data, assuming a spectral index of a = 1.5 (where fu oc v~'^). Column 3 gives the R-band absolute magnitudes (M^j) of the 
nuclear component as derived from the current modelling of the HST data assuming a spectral index of a = 0.2. Columns 4 and 5 list 

the R — K colours of the host galaxy and nuclear component respectively. These colours were derived by combining 12-arcsec aperture 
i?-band photometry from our HST-based models with the 12-arcsec aperture K-haad photometry derived by Taylor et al. (1996) and 
McLure, Dunlop & Kukula (2000), to minimize the uncertainty introduced by errors in constraining the galaxy scale-lengths at K (see 
text for further details). 
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an artifact of the inclusion of RQQ hosts with a substantial 
disc component, since if attention is confined to the 9 RQQs 
with solid elliptical model fits the values remain virtually 
unchanged, i.e.: 

(Mfl) = -23.30 ±0.17 median=-23.30 (RQQ) 

A large variation in the luminosity difference between RQQ 
and RLQ host galaxies has been reported in the literature. 
Differences have ranged from RQQ hosts being fainter than 
their RLQ counterparts by 0.7 — > 1.0 magnitudes in opti- 
cal studies (eg. Smith et al. (1986), Veron-Cetty & Woltjer 
(1990), Bahcall et al. (1997)), to no formal luminosity differ- 
ence being detected from the near-infrared imaging of this 
HST sample (Taylor et al. 1996). It is therefore interesting 
to compare the difference in luminosity detected here with 



the results of the other two recent HST i?-band imaging 
studies of Hooper et al. (1997) and Boyce et al. (1998), both 
of which also used two-dimensional modelling to analyse the 
host galaxies. Below we list the mean differences detected in 
the three studies {{Mrqq) — (Mrlq)) with the associated 
standard error. In calculating the figures for the Hooper et 
al. programme, where no attempt was made to distinguish 
the host galaxy morphologies, the luminosities of the best-fit 
r^^* model have been used. 



AMr = 0.43 ± 0.20 This Work 

AMr = 0.53 ± 0.23 Hooper et al. 

AMr = 0.67 ± 0.29 Boyce et al. 
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Figure 2. Sub-sample histograms of the host-galaxy, absolute, 
integrated il-Cousins magnitudes. 



Figure 3. Sub-sample histograms of the nuclear absolute inte- 
grated ii-Cousins magnitudes. 



The clear implication from these results is that the host 
gala^xies of RQQs are consistently fainter than those of RLQs 
by ~ 0.5 magnitudes in the J?-band. Unlike in earlier stud- 
ies, this difference can no longer be attributed to the model 
fitting of RQQs producing disc fits, since only 2 of the 26 
RQQ host magnitudes included in the above figures are de- 
rived from an exponential host model. The common bias of 
the RLQ redshifts being consistently higher than those of the 
RQQs can also be firmly rejected as a possible cause of the 
difference, with the two quasar types having well-matched 
redshift distributions in all three studies. At least for the 
results presented here, and those of Hooper et al., the selec- 
tion of the RLQ and RQQ quasar samples to have matched 
optical magnitudes, also excludes the possibility that the 
host magnitude difference is as a result of the RLQs being 
intrinsically more luminous. However, it is possible that the 
host-galaxy results of Boyce et al. may have been influenced 
by this effect, considering that the RLQs in that study are 
intrinsically 0.6 ± 0.8 magnitudes brighter than the RQQs, 
although, as can be seen from the large error associated with 
this difference, the overlap of the total quasar luminosities 
is significantly larger than that of the hosts. It is worth not- 



ing at this point that although the results of Taylor et al. 
certainly do not formally support a difference in host mag- 
nitudes (with AMk = 0.4 ± 0.3), they are in fact perfectly 
consistent with the three sets of il-band HST results. The 
question of whether there is any detectable correlation be- 
tween the host-galaxy and nuclear luminosities is explored 
below in Section 6.2.3. 

A general result that can be taken from the luminosities 
of the best-fitting host galaxies is that, with the exception 
of the two Soyfert objects, all of them lie at the extreme end 
of the elliptical galaxy luminosity function. Adopting = 
—22.2 (Lin et al. 1996), after having converted the published 
value of = —21.8 to an integrated magnitude, it can be 
seen that all of the hosts have luminosities of L > L* , with 
25 of the 33 having luminosities L > 2L* . These results 
are in good agreement with those of the previous K-hand 
imaging study which also found all the hosts to be brighter 
than L^. Independent support for this result comes from 
the findings of Hooper et al. (1997) and Boyce et al. (1998), 
both of which used the F702W (wide R) filter on WFPC2. 
Using the Lin et al. value for M^, Hooper et al. found 15 of 
their 16 {z ~ 0.4) quasars to have L > L*, while Boyce et 
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al. found all 11 of their 14, objects for which a host model 
was fitted, to have L > L* (both sets of results having been 
converted to our adopted cosmology). 



6.2.2 Nuclear Luminosities 

The luminosities of the best-fitting nuclear components are 
listed in Tabic 5 and illustrated in Fig 3, with sub-sample av- 
erages summarized below. The radio galgixy 0230—027 has 
been excluded from these summary statistics as the mod- 
elling code detected no unresolved component in this object. 



(Mr) 
(Mr) 
{Mr) 
(Mr) 



-22.27 ±0.45 
-19.27 ±0.64 
-24.07 ±0.22 
-22.95 ±0.57 



median=— 23.36 
median= — 19.64 
median=— 24.01 
mcdian=-23.70 



(ALL) 
(RG) 
(RLQ) 
(RQQ) 



As expected, the unresolved nuclear components displayed 
by the RGs are nearly two orders of magnitude weaker than 
those of the RLQs, in good agreement with the orientation- 
based radio-loud unification scheme first described by Pea- 
cock (1987) and Barthel (1989) (see Section 6.7). It would 
appear from the mean figures given above that there is a 
substantial difference in the nuclear components of the RQQ 
and RLQ quasar samples. However, the noticeable offset in 
the difference between the mean and median figures, com- 
bined with cxaminatiou of histogram shown in Fig 3, reveals 
that the RQQ mean is being biased by the three objects 
which have total luminosities fainter than My = —23.0. If 
these throe objects arc excluded from the RQQ sample, the 
mean and median values of nuclear luminosity become: 



(Mr) 



-23.95 ± 0.26 median=-23.97 (RQQ) 



completely consistent with the equivalent RLQ figures. The 
similarity between the nuclear components of the RQQ and 
RLQ sub-samples is reassuring considering that the sample- 
selection process was indeed originally designed to produce 
two quasar samples well-matched in luminosity. This result 
also re-affirms that the measured difference in the respec- 
tive host magnitudes of the RQQs and RLQs cannot be at- 
tributed to some bias arising from the RLQs in our sample 
having, on average, more powerful nuclei. 



6.2.3 Quasar Host-Nuclear Luminosity Correlation 

Because the measured difference between the RQQ and RLQ 
host-galaxy magnitudes is substantially greater than the 
average difference between their respective nuclear compo- 
nents (excluding the three sub-luminous objects) it is not 
expected that there should be a strong host-nuclear lumi- 
nosity correlation. This is indeed the case, aud application 
of the Spearman Rank correlation test confirms that there 
is no evidence for any correlation, returning a probability 
p = 0.95 that the null hypothesis of no correlation is accept- 
able. 

The null result found here contrasts with the positive 
correlations found previously by many authors (eg. Smith 
et al. (1986), Bahcall et al. (1997), Hooper et al. (1997)). 
The suspicion often cast upon the reality of a positive liost- 
gala^xy: nuclear correlation is that there are two obvious se- 
lection effects at work. The first of these is that it is ob- 
viously much more difficult to detect faint galaxies which 



host brighter quasars. It seems clear that given the difficul- 
ties associated with ground-based seeing, and the perils of 
PSF-subtraction on HST data (eg. Bahcall et al. 1997), that 
this could well be a contributing factor. In the present study 
there are several RQQs with bright nuclear components, and 
relatively faint host galaxies (eg. 0953±415), providing con- 
fidence that the techniques employed here have been succes- 
ful in overcoming this possible source of a false correlation. 
The other selection effect which could contribute to a false 
positive result is that weak AGN in bright host galaxies will 
not be classified as quasars. 

Although the results presented here for powerful \ow-z 
quasars show no evidence for a correlation, it docs appear 
to exist for lower-luminosity Seyfert galaxies. In the case 
of Seyferts, the reduction in both redshift and contrast be- 
tween the host and nuclear components means that it can 
be confidently assumed that all of the host galaxies are be- 
ing detected, removing one of the sources of bias mentioned 
above. Indeed, the work of McLeod & Rieke (1994) on two 
samples of low- and high-power AGN at low-2; has led them 
to suggest that the relation between host and nuclear lu- 
minosity is of the form that there is a minimum host lu- 
minosity required to produce a particular quasar luminos- 
ity. This would be consistent with both the positive cor- 
relation found at lower AGN power, and the fiat relation 
found here for higher powered quasars. The host galaxies 
for the quasars studied here have already been shown to be 
among the brightest known, and therefore demand that the 
galaxy-quasar relation tails-off at high quasar luminosity. 
The question of how the nuclear and host-galaxy luminosi- 
ties relate to central black-hole mass and quasar fuelling rate 
is pursued in Section 8. 

6.3 Scale-lengths 

The best-fitting values of the host-galaxy scale-lengths listed 
in Table 3 are displayed as histograms for the separate sub- 
samples in Fig 4. The average figures for the three sub- 
samples are : 



(ri/a) = 12.23 ± 1.00 
{t-i/2> = 13.76 ± 2.18 
(ri/a) = 11.73 ± 1.07 
{ri/2> = 11-45 ± 1.66 



median=10.45 
median=11.27 

median=11.28 
median= 9.30 



(ALL) 
(RG) 
(RLQ) 
(RQQ) 



Taking the median figures, the RG and RLQ sub-samples 
can again be seen to be in remarkable agreement, while the 
RQQ hosts are ~ 20% smaller. 

Galaxy half-light radii of ~ 10 kpc mark out these radio 
galaxies and quasar hosts as being giant elliptical galaxies. In 
a study of the galaxies in the Virgo cluster, Capaccioli et al. 
(1992) found that beyond a scale- length of ~ 3 kpc the only 
galaxies to be found were extremely luminous ellipticals, or 
cD type brightest cluster galaxies (BCG). The question of 
how the properties of the host galaxies studied here relate 
to those of BCGs is explored further in Section 7. In further 
support of the radio-loud unification scheme the scale-length 
figures obtained here for the RG and RLQ sub-samples can 
be compared with the B- and F-band study of Smith & 
Hcckman (1990), who determined the scale- lengths of 41 
powerful radio gala^xies in the redshift range < z < 0.26. 
The median scale-length figure determined is equivalent to 
~ 17 kpc in the cosmology adopted here, somewhat larger 
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Figure 4. Histograms of the best- fitting host-galaxy scale-lengths 
for the three AGN sub-samples. 

than the results presented above. However, if the comparison 
is restricted to the 22 objects studied by Smith & Heckman 
which have strong optical emission lines, more likely to have 
the same FR II morphology as the majority of radio-loud 
objects studied here, then the median scale-length falls to 
13 kpc, in excellent agreement with the figures presented 
above. 

If the host scale-lengths arc plotted against the inte- 
grated model luminosities, as in Fig 5, a tight correlation 
is found, and a least-squares fit to these data produces a 
relation of the form L oc r'^''^'' . This is in impressively close 
agreement with the corresponding result L ocr'^'^^ found for 
I0W-2: inactive ellipticals by Kormendy (1977). Given this 
result, it is obviously of interest to ask whether the best-fit 
model parameters for /ii/2 and ri/2 given in Table 3 do in 
fact yield a Kormendy relation of slope ~ 3. 

6.4 The Kormendy Relation 

The /ii/2 — ri/2 relation for the host galaxies of all 33 objects 
is shown in Fig 6. For the four RQQ sources which have 
combined disc/bulge fits it is the parameters for the bulge 
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Figure 5. Plot of absolute host magnitude against best- fitting 
host scale- length for the RGs (crosses), RLQs (open circles) and 
RQQs (filled circles) in the HST sample. The solid line is the least- 
squares fit to the data, which has the form L (x r'''^^ . For the four 
RQQs which required a combined disc/bulge fit, the best-fitting 
bulge parameters have been plotted. 



component that have been plotted. The least-squares fit to 
the data (solid line) has the form: 



Ml/2 = 2.90±o.22 logio ^1/2 + 18.35dbo.22 



(ALL) 



(errors arc ±lcr) showing for the first time from optical imag- 
ing, that the host galaxies of quasars lie on this photometric 
projection of the fundamental plane. The individual fits to 
the /11/2 — ri/2 relations for the individual sub-samples are: 

fly 2 = 2.86±o.32 log ri/2 + 18.44±o.36 (RG) 
^1/2 = 3.98±o.7i log ri/2 -I- 17.02±o.7B (RLQ) 
Ml/2 = 2.99±o.34 log ri/2 -I- 18.39±o.30 (RQQ) 

It is readily apparent from the individual Kormendy re- 
lations that the RG and RQQ sub-samples follow very simi- 
lar relations, consistent with each other in terms of slope and 
normalization. One encouraging aspect of this is that the 
best-fitting bulge components for the RQQ objects which 
have combined disc/bulge model fits lie naturally on the 
Kormendy relation defined by the other nine, single compo- 
nent, RQQ hosts. With no restrictions placed on the range of 
parameter values available to the modelling code, the fact 
that the fitted bulge components have physically sensible 
values gives further confidence that the introduction of the 
combined disc/bulge fits was justified, and does not involve 
over-fitting of the data. 

It would appear from the individual fits given above 
that the RLQ sub-sample follows a steeper slope, compared 
with the best-fit relation to the other two sub-samples. It 
is also noticeable that the 1 a error returned by the least- 
squares fitting procedure is more than twice that returned 
from the fitting to the RG sub-sample (which also has 10 
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Figure 6. The Kormendy relation followed by the hosts of all 

33-objccts in the HST sample. For the four RQQs which have a 
significant disc component the best-fitting bulge component has 
been plotted. The solid line is the least-squares fit to the data 
which has a slope of 2.90, in excellent agreement with the slope 
of 2.95 found by Kormendy (1977) for inactive ellipticals in the 
S-band. The dotted line has a slope of 5, indicative of what would 
be expected if the scale-lengths of the host galaxies had not been 
properly constrained (see text) 



objects), suggesting that there may well be outliers within 
the RLQ sub-sample which axe biasing the fit. 

This is indeed the case, with two outlying RLQ hosts 
(1004+130, 2141+175), one of which (1004+130) is the most 
nuclear-dominated RLQ in the sample, apparently distort- 
ing the fit due to the probable under-estimation of their 
host-galaxy scale-lengths. If these two objects are removed 
the best-fitting relation for the RLQ hosts becomes: 



/il/2 = 3.19±0.67log 



10 '1/2 



+ 17.95±o.72 



(1) 



Thus, with the exclusion of the two outlying RLQ hosts, 
the three Kormendy relations followed by the sub-samples 
are internally consistent with each other and, perhaps more 
importantly, formally inconsistent with a slope of 5 which, 
as highlighted by Abraham et al. (1992), would be a clear 
indication that the data and/or the modelling procedure 
employed to analyse the images is unable to break the /^l/2 — 
ri/2 degeneracy. 

We stress this point here because some authors have 
claimed to have demonstrated that AGN hosts follow a Ko- 
rmendy relation, when in fact the derived relation simply 
reflects a failure to constrain galaxy scale-lengths. Indeed 
some authors (e.g. de Vries et al. 2000) have even mistak- 
enly interpreted a slope of 5 as evidence that AGN hosts 
do not follow the same Kormendy relation as normal giant 
ellipticals. 

The reason that a 'pseudo Kormendy relation' with a 
slope of 5 can be produced with inadequate data and/or 



modelling follows simply from the fact that the integrated 
luminosity of a galaxy model follows the relation: 



OC Il/2ri/2 



(2) 



where I1/2 is the surface-brightness at ri/2- Given that 
Ml/2 OC — 2.51og7i/2, it follows that if the modelling pro- 
cedure can successfully constrain the host luminosity, but 
not the scale-length, then the apparent best-fitting values of 
/ii/2 and ri/2 will be randomly distributed along a relation 
obeying: 



/^l/2 OC 51ogri/2 



(3) 



with appropriate normalization to fit the integrated luminos- 
ity. Due to the fact that the host galaxies display a relatively 
small spread in luminosity (mean Mr = —23.44, cr = 0.64) 
it might be expected that, if our analysis was unable to 
accurately determine the individual host scale-lengths, the 
resulting /ii/2 — ''1/2 relation would be well fitted with a slope 
of 5 (as found by de Vries et al. 2000) and a normalization 
to match = —23.44. Such a relation has the form: 



Ml/2 = 5.0 logio ^1/2 + 16.40 



(4) 



and is plotted as the dotted line in Fig 6. It is obvious from 

Fig 6 that this relation is not consistent with the data, giving 
confidence that the methods of analysis employed here have 
allowed the accurate determination of the host-galaxy scale- 
lengths. 

The diagnostic power of Fig 6 is discussed in the context 
of recent studies of the fundamental plane in Section 7.2. 



6.5 Axial Ratios 

If the AGN host galaxies are indeed indistinguishable from 
massive inactive ellipticals then they should display an ax- 
ial ratio distribution which is also identical. The axial ratio 
distribution of normal ellipticals is well studied and known 
to peak at values of h/a > 0.8 (Sandage, Freeman & Stokes 
1970; Ryden 1992). The best-fitting axial ratios from the 
two-dimensional modelling are listed in Table 3 and dis- 
played in the form of separate sub-sample histograms in Fig 
7. The corresponding mean and median values are: 



{b/a) 


= 0.81 + 0.02 


median^ 


=0.85 


(ALL) 


{h/a) 


= 0.86 + 0.03 


median^ 


=0.90 


(RG) 


{b/a) 


= 0.80 ± 0.05 


median= 


=0.82 


(RLQ) 


{h/a) 


= 0.78 ± 0.03 


median= 


=0.84 


(RQQ) 



It can be seen from these figures that the axial ratio distri- 
butions of all three sub-samples are perfectly consistent with 
that expected from a sample of elliptical galaxies. There is a 
slight suggestion that the RGs have higher axial ratios than 
average, but this is not formally significant. 

The axial ratio results from the two-dimensional mod- 
elling agree well with the recently-published findings of 
Boyce et al. (1998) ,who found that all 11 of the quasars 
from their 14-object (z ~ 0.3) sample for which a model fit 
was possible, displayed axial ratios with h/a > 0.65. In con- 
trast. Hooper et al. (1997) found only 2 of their 16 {z ~ 0.45) 
quasars to have axial ratios with b/a > 0.6. The reasons for 
this apparently contradictory result probably lie in a com- 
bination of the higher redshifts of the Hooper et al. objects, 
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Figure 7. The axial ratio distributions for the three host-galaxy 
sub-samples as determined by the two-dimensional modelling. 

and their use of the PC instead of the WF detectors. As 
has been explained, these two factors will undoubtedly re- 
sult in a reduced sensitivity to low surface-brightness fear 
tures, given that the exposure times used were the same 
as adopted in this study. In fact, Hooper et al. noted their 
concern that their modelling technique of two-dimensional 
cross-correlation could have been influenced by high surface- 
brightness features such as bars or tidal tails. Given the clear 
result presented here, and the support of the Boyce et al. 
results at similar redshifts, this seems the most likely expla- 
nation of the Hooper et al. result. 

6.6 Host-Galaxy Colours 

The results presented in the last five sub-sections strongly 
suggest that in terms of morphology, luminosity, scale- 
length, Kormendy relation, and axial-ratio distribution the 
hosts of powerful AGN are identical to normal, massive, in- 
active ellipticals. The one final parameter which can be read- 
ily recovered from the modelling is host-galaxy colour. The 
desire to obtain reliable optical-infrared colours for the host 
galaxies was one of the original motivations for this HST 



o 

12 3 4 

R-K 

Figure 8. The rest- frame R—K colours for the three host-galaxy 
sub-samples. The three sub-samples can be seen to be consistent 
with each other, tightly distributed around a value of R—K ~ 2.5. 

imaging study. Given that the results presented thus far sug- 
gest the host galaxies are otherwise normal massive ellipti- 
cals, it is now even more interesting to investigate whether 
the hosts also display the red colours associated with old 
stellar populations, or whether they have significantly bluer 
colours, indicative of either a generally young stellar pop- 
ulation, or of substantial secondary star formation induced 
by interactions, the central AGN, or both. 

The simplest and most obvious way to calculate R ~ K 
colours for the host galaxies is to combine the integrated 
optical magnitudes presented in Table 5 with the integrated 
absolute /f-band magnitudes derived for the sample by Tay- 
lor et al. (1996). However, there are two potentially serious 
problems with this strategy. First, Taylor et al. did not al- 
ways decide on the same host morphology as has now been 
revealed by the HST imaging, and so in some cases their 
quoted absolute K magnitude is derived from, for example, 
a disc fit which originally appeared marginally preferable 
to a de Vaucouleurs law. This obviously matters because 
the disc fits to the JsT-band data were 0.5 — » 1.0 magni- 
tudes fainter than the equivalent elliptical fits. Second, in 
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the light of the new, better constrained models derived from 
the HST data, and the results of more recent tip-tilt infra- 
red imaging reported by McLure, Dunlop & Kukula (2000), 
it is now clear that the half-light radii derived by Taylor et 
al. wore systematically over-estimated (most likely duo to 
the '--J 1" seeing and coarse spatial resolution (0.62"/pix) of 
the original K-band observations - see also Simpson et al. 
2000). Given that the integrated luminosity of the standard 
de Vaucouleurs r^^* law is proportional to /i/2'"ei it is clear 
that scale-length errors of the order of 1.5 — 2 could seriously 
bias the integrated K-h&nd luminosity. 

Therefore, to obtain a robust value for the R—K colour 
of each host from the existing data we proceeded as follows. 
First, the K-band luminosities were based on the best-fitting 
model with the same morphology as the equivalent HST re- 
sult. Second, the R — K colours were based on 12"-diameter 
aperture photometry performed on both the HST and K- 
band model fits, to minimize the impact of scale-length er- 
rors on the derived colour, while at the same time including 
the vast majority (> 75%) of host-galaxy light. 

The results of these host-galaxy R — K calculations 
are listed in Table 5, with absolute R — K histograms for 
the separate sub-samples shown in Fig 8. The absolute 
colours have been calculated assuming spectral indices of 
a = 1.5 & a = 0.0 for the R- and if-bands respectively. 

It is readily apparent from Fig 8 that the colour dis- 
tributions of the three sub-samples are consistent with each 
other, with mean and median absolute colours: 
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Figure 9. Plot of apparent R — K colour versus absolute i?-band 

luminosity of the fitted nuclear components of the objects in the 
HST sample. Shown in the plot are the RG nuclei (crosses), RLQ 
nuclei (open circles) and RQQ nuclei (filled circles). Also shown 
is the least-squares fit to the data which has a slope of 0.68±0.05. 



(R-K)= 2.48 ± 0.05 
(R-K) =2.47 ±0.09 

{R- K) = 2.48 ±0.07 
{R- K) = 2.48 ± 0.08 



median=2.44 
median=2.56 

mcdian=2.46 
median=2.41 



(All) 
(RG) 

(RLQ) 
(RQQ) 



As demonstrated by Nolan et al. (2001) such colours are con- 
sistent with those expected of an evolved stellar population 
of age 10 — 13 Gyr. The homogeneity of these colours pro- 
vides a further strong indication that the hosts of all throe 
classes of AGN are derived from the same parent population 
of massive, well-evolved elliptical galaxies. The inevitable 
conclusion from this is that these galaxies, or at least the 
bulk of their stellar populations, must have formed at high 
redshift, a theme which is revisited in Section 9. The fact 
that the R — K colours are so similar to passive ellipticals 
forces the further conclusion that any star-formation asso- 
ciated with AGN activity must either be dust enshrouded, 
confined to tidal features masked from the modelling, or pos- 
sibly restricted to the central ^ 1 kpc of the host (in which 
case the additional emission would almost certainly be at- 
tributed to the unresolved nuclear component). Large-scale 
star-formation involving a substantial fraction of the mass 
of the host galaxy, is effectively ruled out by these results. 



6.7 Nuclear Colours: Implications for Unified 
Models of Radio-Loud AGN 

The R — K colours of the fitted rmclei for all 33 objects (ex- 
cept the RG 0230—027 - see above) are plotted against abso- 
lute nuclear R magnitude in Fig 9. It is clear from this plot 
that there is no significant difference between the colours of 



the nuclei of the RQQs or RLQs. More importantly, how- 
ever, this plot allows us to explore the relationship between 
the relatively weak nuclear components found in the RGs, 
and the (naturally) brighter components found in the RLQs. 

This is of value because it allows us to perform another 
test of the proposed unification of RLQs and RGs by ori- 
entation (Peacock 1987, Barthel 1989). As discussed above, 
our host-galax;y results can already be viewed as providing 
strong support for such a picture, due to the fact that the 
basic properties (i.e. luminosity, size, axial ratio, age etc) of 
the RGs and RLQ hosts in our matched samples are indis- 
tinguishable. However, it must also bo realised that identical 
RG and RLQ host galaxy properties would also be expected 
in a scenario where the two classes of radio-loud AGN are in 
fact linked by evolution in time, rather than by orientation. 

Fig 9 offers the possibility of differentiating between 
these two alternative explanations, by allowing us to test 
whether the R—K colours of the RG nuclear components are 
consistent with dust-obscured quasar nuclei. If it is assumed 
that the reddening of the RG nuclei obeys a k,\ cc law, 
then the reddening vector which should connect the RGs and 
RLQs in Fig 9 would be expected to be given by E{R—K) = 
0.70Afl. 

In fact the slope of the least-squares fit in Fig 9 (solid 
line) is 0.68 ± 0.05, which corresponds to a dust reddening 
law kaA""'®^^"'^^. Furthermore, if the fitted _R-band nuclear 
components of the RGs are de-reddened using a slope of 
0.68, until each object has the mean colour of the RLQ nuclei 
{R — K = 3.0), then the mean luminosity of the de- reddened 
RG nuclei is Mr, = -23.90 ± 0.46. It can be seen that this 
is perfectly consistent with the previously determined mean 
nuclear luminosity of the RLQ sample, Mn = — 23.73±0.10. 
Given the clean nature of this result, there seems little room 
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for argument that the RG and RLQ sub-samples are most 
consistent with being drawn from the same population of 
objects, viewed from different orientations, with every RG 
displaying evidence of a buried quasar nucleus. 



7 THE RELATION OF QUASAR HOSTS TO 
'NORMAL' GALAXIES 

7.1 Dependence of Host Morphology on Nuclear 
Luminosity 

In the preceding section wc presented a scries of results 
which together lead to the inescapable conclusion that the 
hosts of virtually all powerful AGN are essentially normal, 
massive ellipticals. This result is particularly clean for the 
radio-loud objects in our sample. However, exceptions to any 
rule are often of great significance, and one should not ignore 
the fact that we have found two disc-dominated hosts in our 
RQQ sample, along with two bulge-dominated RQQ hosts 
which do, however, possess significant disc components. 

What might such exceptions to the general 'perfect el- 
liptical' rule be telling us? An important clue comes from the 
fact that, as discussed above, the two disc-dominated hosts 
in our sample transpire to contain the two least-luminous 
nuclei. Indeed, in much of the preceding analysis, these two 
objects have had to be deliberately excluded from sample 
comparisons due to the fact that they are not actually lumi- 
nous enough to be classed as quasars. However, the fact that 
these two radio-quiet low-luminosity interlopers are the only 
RQQs in our sample with disc-dominated hosts strongly sug- 
gests that host morphology is a function of nuclear optical 
luminosity. 

Our own samples do not span a sufficient range in op- 
tical luminosity to allow a statistical test of this hypothesis. 

However, we can explore this possibility, and place our re- 
sults in a wider context, by combining them with those of 
Schade et al. (2000) who have studied the morphologies of 
the host galaxies of a large sample of X-ray selected AGN 
spanning a wider but lower range of optical luminosities. 
This we have done in Fig 10, which shows the ratio of bulgc- 
to-total host luminosity plotted as a function of nuclear op- 
tical power. 

This instructive diagram illustrates a number of impor- 
tant points. First, it shows that while pure spheroidal galax- 
ies can host AGN with optical luminosities ranging over 
two orders of magnitude, very disc-dominated hosts start to 
die out for Mv{nuc) < —21 and no disc-dominated galax- 
ies appear capable of hosting nuclear emission more lumi- 
nous than Mv — —23 (close to the traditional, albeit rather 
ad hoc, quasariseyfert classification boundary). Second, this 
plot confirms that the two most 'discy' hosts uncovered in 
the present study do indeed lie in a region of the parameter 
space where disc-dominated host galaxies are quite common. 
Third, in the context of this diagram, the persistence of sig- 
nificant discs around two of the more luminous RQQs in our 
sample (0052H-251 and 0157-t-OOl) seems perfectly natural, 
and certainly not inconsistent with the clear trend towards 
more spheroid-dominated hosts with increasing nuclear out- 
put. 

Fig 10 clarifies the origin of much of the confusion which 
has surrounded the nature of RQQ hosts. In particular it 
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Figure 10. Plot of relative contribution of the bulge/spheroidal 
component to the total luminosity of the host galaxy versus ab- 
solute V-band luminosity of the fitted nuclear components. The 
plot shows the results for our own HST sample (RLQs as open 
circles, RQQs as filled circles) along with the results from Schade 
et al. (2000) for a larger sample of X-ray selected AGN (stars). 
This plot illustrates very clearly how disc-dominated host galax- 
ies become increasingly rare with increasing nuclear power, as is 
expected if more luminous AGN are powered by more massive 
black holes which in turn are housed in more massive spheroids. 



is clear that studies of RQQ samples which extend to in- 
clude significant numbers of objects with nuclei less lumi- 
nous than Mv ^ —23 are likely to uncover significant num- 
bers of hosts with substantial discs, in contrast to the results 
of the present study. This does indeed appear to be the case. 
For example, Bahcall et al. (1997) and Hamilton et al. (2002) 
have reported that approximately one third to one half of 
radio-quiet quasars lie in disc-dominated hosts. However, if 
one restricts their results to quasars with nuclear magni- 
tudes Mv < -23.5 (for which wc find 10 out of 11 RQQs 
lie in ellipticals), it transpires that 6 out of the 7 remaining 
objects in the Bahcall et al. study lie in ellipticals, while 
17 out of the remaining 20 objects in the archival study of 
Hamilton ct al. also have elliptical host galaxies. 

Fig 10 is also at least qualitatively as expected if AGN 
are in fact drawn at random from the general galaxy popula- 
tion, subject to the constraint that the host galaxy contains 
a black hole of sufficient mass to produce the observed nu- 
clear output (see also Wisotzki, Kuhlbrodt & Jahnke 2001). 
As discussed further in Section 8, a black-hole mass nibh > 
5 X 10* Mq appears to be required to produce a quasar with 
Mv < —23.5. When this constraint is combined with the 
now apparently inescapable result that black-hole mass is 
proportional to spheroid mass Mih — 0.001 — 0.003Mspft, 
then Fig 10 can be viewed as a simple manifestation of the 
fact that spheroids in the present-day universe with baryonic 
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masses > 3 x 10^ ^Mq are rarely accompanied by significant 
discs (and are thus classed as massive ellipticals). 

This conclusion leads naturally to the question of why 
most massive black holes in the present-day universe should 
be inactive, while a subset arc receiving sufficient fuel to 
shine as quasars. Below we present the results of a series 
of statistical comparisons designed to decide whether there 
are in fact any significant observable differences between the 
hosts (and host environments) of quasars and comparably 
massive, but inactive, elliptical galaxies. 

7.2 Comparison with ULIRG merger remnants 

Before proceeding to compare quasar hosts with old 
passively-evolving massive ellipticals, it is worth pausing 
to consider whether host galaxies with the properties we 
observe could plausibly be produced as the outcome of a 
(relatively recent) violent merger of two disc galaxies. This 
might already seem fairly implausible given the series of 
results presented in Section 6, in particular the red host- 
galaxy colours consistent with old evolved populations. Nev- 
ertheless, it is well known from both simulations (Barnes 
& Hernquist 1992) and observations that the red^infrared 
light from ULIRGs often follows an r'^'^* law (Wright et al. 
1990, Scoville et al. 2000) and it is possible that dust extinc- 
tion in such objects could also produce red colours (albeit 
such colours would be unlikely to mimic so well the proper- 
ties of an evolved population, as discussed by Nolan et al. 
2000). Moreover, several authors have proposed an evolu- 
tionary scheme in which ULIRGs might evolve into radio- 
quiet quasars on a time-scale < 1 Gyr (Sanders & Mirabel 
1996). 

Fortunately, the data on quasar hosts presented here 
are of sufficiently high quality to completely exclude the 
possibility that they are the recent outcome of the violent 
disc-disc mergers which appear to be the origin of most 
ULIRGs in the low-redshift universe. While it is true that 
ULIRGs such as Arp220 have surface brightness profiles 
well-described by an r^/*-law, recent observations (Genzel 
et al. 2001) have shown that such remnants lie in a com- 
pletely different region of the fundamental plane than that 
which we have shown is occupied by quasar hosts. Specifi- 
cally, the effective radii of ULIRGs is typically ~ 1 kpc, an 
order of magnitude smaller than the typical effective radii 
of the quasar hosts as illustrated in Figs 4, 5 and 6. 

Indeed, one can go further and conclude that whereas 
ULIRGs can be the progenitors of the population of 
intermediate-mass ellipticals which display compact cores 
and cusps (Faber et al. 1997), the quasar hosts lie in a re- 
gion of the /ie — re plane which is occupied by boxy, giant 
ellipticals with large cores, a large fraction of which lie in the 
centres of clusters (Faber et al. 1997). This is demonstrated 
in Fig 11 where we have constructed a composite fl-band 
Kormendy diagram, combining our own results on the 31 
bulge-dominated AGN hosts with the data from Genzel et 
al. (2001) on LIRGs and ULIRGs, and the data from Faber 
et al. (1997) on 'discy' and 'boxy' ellipticals. 

This plot provides strong evidence that the quasar hosts 
belong to the class of largo 'boxy' ellipticals, which vari- 
ous authors have suggested form at earlier times than their 
lower-mass 'discy' counterparts (Kormendy & Bender 1996, 
Faber et al. 1997). 
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Figure 11. A comparison of the properties of the AGN hosts 
with those displayed by various other types of spheroid on the 

photometric projection of tlic fundamental plane. Symbols for 
the quasar hosts and radio galaxies are as in Figs 5 and 6. The 
stars arc the data for ULIRGs and LIRGs from Genzel et al. 
(2001) transformed from the infrared to the _R-band assuming 
R — K=2.5. Triangles and squares indicate the positions of 'discy' 
and 'boxy' ellipticals from Faber et al. (1997) after conversion to 
Ho = 50 kms-iMpc-i. 

In summary, all the available evidence indicates that 
luminous low-redshift quasars are the result of the re- 
triggering of a massive black hole at the centre of an old 
evolved elliptical, and do not generally occur in the rem- 
nants of the recent disc-disc mergers which power ULIRGs. 
While there is one RQQ in our sample which is also a ULIRG 
(0157-1-001), in this case it is clear that at least one of the 
merging galaxies already has a massive spheroid capable of 
containing a pre-existing massive black hole. Thus, while 
there is clearly some overlap between the ULIRG and quasar 
phenomenon, present evidence suggests that, at least at low 
redshift, the ULIRG — > quasar evolutionary route can only 
apply to a fairly small subset of objects. 

7.3 Comparison with Brightest Cluster Galaxies 

Having established that the hosts of quasars are massive 
ellipticals, it is then of interest to explore how their proper- 
ties compare to those of inactive comparably-massive ellip- 
ticals at similar redshift. In this sub-section we compare the 
properties of the quasar hosts, and in particular their scale- 
lengths, with those of bright galaxies found in rich clusters. 
This comparison is motivated, in part, by the location of the 
host galaxies on the /Xe — re projection of the fundamental 
plane as discussed in the previous subsection. However, this 
motivation is further strengthened by the results of a study 
of the environments of the quasars and radio galaxies un- 
dertaken by McLure & Dunlop (2000). McLure & Dunlop 
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(2000) used the full WFPC2 images of our AGN sample to 
measure the spatial clustering amplitude (Bgg). This quan- 
tity is inevitably rather poorly constrained for each individ- 
ual object, but, on average, they found that all three classes 
of AGN typically inhabit environments as rich as AbcU Class 
~ 0. Moreover, several objects appeared to lie in noticeably 
richer environments (Abell Class 1). 

Given these clustering results, and the fact that the 
quasar hosts are clearly comparable to many BCGs in terms 
of optical luminosity, it is interesting to explore whether 
the sizes of the quEisar hosts are as expected in the light of 
their apparent cluster environment. Strong circumstantial 
evidence that the quasar host galaxies arc directly compara- 
ble to first-ranked cluster galaxies comes from a comparison 
of the host-galaxy Kormendy relation with that found for 
cD galaxies by Hamabc & Kormendy (1987) in the B-band: 

/.i/2 = 2.94 logiori/2 + 20.75 (5) 

If a typical elliptical gala^Ky colour a.t z = 0.2 oi B — R = 2.4 
is assumed (Fukugita et al. 1995), then the inferred B-band 
Kormendy relation for all 33 objects in the HST sample 
presented in Section 6.4 becomes: 

/xi/2 = 2.90 logio '"1/2 + 20.75 (6) 

which can be seen to be in excellent agreement with the 
Hamabc & Kormendy result. 

To explore further the connection between host clus- 
ter environment and galaxy size, a sample of 51 first-ranked 
cluster galaxies (confined to the same redshift range as the 
HST sample but spanning a range of Abell classes) was con- 
structed from the study of Hoessel & Schneider (1985). Ap- 
plication of the KS test to the scale-length distributions of 
the AGN host galaxies, and the 14 first-ranked cluster galax- 
ies in this sample from Abell clusters of class & 1, shows 
that the two are not significantly different {p = 0.22). How- 
ever, the extension of the comparison sample to include the 
first-ranked gala^xies of Abell class 2 clusters (15 objects), 
shows these two distributions to be different at the 2cr level 
(p = 0.011). 

This is an very interesting result because it provides 
independent support for the results of the cluster analysis of 
McLure & Dunlop (2000) which also favour an environment 
more consistent with Abell class 0, or at most (in only a few 
cases) Abell class 1. 

Finally, to check that the ground-based results discussed 
above can be confirmed by modelling of more recent HST 
data, and to minimize possible systematic errors, we decided 
to model a series of HST images of bright cluster galaxies 
(at comparable redshift to our quasars) using the same 2-D 
modelling code used to model the AGN hosts. HST WFPC2 
images of five Abell clusters with similar redshifts to the ob- 
jects in the HST sample were therefore retrieved from the 
HST archive facility^ . Following the production of reduced 
images of comparable signal-to-noise to the host galaxy im- 
ages, the four brightest galaxies in each cluster were iden- 
tified and modelled in identical fashion to the AGN hosts. 
The results of the elliptical model fits to these 20 galaxies 
are presented in Table 6. 

The scale-lengths and absolute magnitudes listed for the 

^ http: / / archive.stsci.edu 



BCGs in Table 6 can be instructively compared with the 
corresponding quantities for the AGN hosts listed in Tables 
3 and 5. First, it is clear that none of our AGN hosts is 
comparable in size to either the first or second ranked clus- 
ter galaxies in very rich, Abell Class 4 clusters. However, 
the largest of the AGN hosts are comparable in size to the 
1st or 2nd ranked galaxies in the Abell Class and Class 
1 clusters. This more detailed galaxy-by-gaJaxy comparison 
thus provides further support for the basic conclusions of 
McLure & Dunlop (2000) already mentioned above. Indeed, 
at least within the radio-loud subsamples, there is a sugges- 
tion of a rather clean link between the biggest galaxies and 
the richest cluster environments. Specifically, the two ra- 
dio galaxies found by McLure & Dunlop (2000) to have the 
richest environments (1.342-016 (see Fig 12) & 0917+459) 
transpire (see Table 3) to be two out of the three RGs with 
re > 20 kpc. Moreover, both the values of Bgg for these ob- 
jects (~1500 and ~1000) and their half-light radii both point 
towards the same conclusion (given Table 6), namely that 
these most massive hosts are the first- or second- ranked 
BCGs in clusters of richness Abell Class 0— »1. A similar 
remarkably clean correspondence can be found within the 
RLQ sample. The two RLQs found by McLure & Dun- 
lop (2000) to have the richest environments (2349—014 and 
0137-1-012) transpire (see Table 3) to be the only two RLQs 
with hosts for which re > 14 kpc. However, no such clean 
correspondence is evident within the RQQ sample. This may 
mean that RQQs, even when found in cluster environments, 
are less likely to be central cluster galaxies than their radio- 
loud counterparts. However, given the large errors on Bgq 
such a conclusion may perhaps be premature. 

In summary, this comparison of the properties of the 
hosts of powerful AGN and bright cluster galaxies supports 
the basic conclusion reached by McLure & Dunlop (2000) 
via the completely independent measure of clustering am- 
plitude. The larger AGN hosts have properties comparable 
to first or second ranked BCGs in Abell class 0^1 clusters. 
Both radio-loud and radio-quiet quasars thus appear to oc- 
cur in a range of environments which is certainly consistent 
with the range of environments occupied by inactive ellip- 
ticals of comparable size and luminosity. While none of our 
host galaxies are comparable to the largest known BCGs, we 
must caution that this does not mean AGN avoid such en- 
vironments (as has sometimes been inferred for FRII radio 
sources from the results of, for example, Prestage & Peacock 
1988). In fact, because Abell Class 3 and 4 clusters are very 
rare, with a number density less than < 1.6 x 10"'' h'^Mpc"^ 
(Croft et al. 1997), even an all-sky survey covering the red- 
shift band of this study (0.1 < z < 0.25) would be expected 
to contain < 1500 such rich clusters. Consequently, given 
that only 1 in ~ 1000 — 10000 massive black holes appears 
to be active in the present-day universe (see Section 9 and 
Wisotzki, Kuhlbrodt & Jahnke 2001), we should not be sur- 
prised by the fact that no low-redshift quasar has yet been 
found to lie within such a rich environment or hosted by an 
ultra-massive BCG. 

7.4 Interactions 

It is often stated in the literature (eg. Smith et al. 1986, 
Hutchings & Neflf 1992, Bahcall et al. 1997) that morpho- 
logical disturbance is a common feature of the host galaxies 
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Table 6. The model results for the four brightest galaxies in each of five clusters obtained from the HST archive, spanning a range in 
richness from Abell Class to Abell Class 4. The redshifts of the clusters are 0.11, 0.19, 0.23, 0.175 & 0.171 respectively. The images 
of A0018, A0103 and A2390 were taken through the F702W (wide R) filter. The images of A2390 and A1689 were taken through the 
F814 (/-band) filter. The conversion between /- and i?-magnitudes has been performed assuming R — I = 0.8 (Pukugita et al). The final 
reduced images had signal-to-noise comparable to the HST host-galaxy data. 
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Figure 12. The full WF2 image of the radio galaxy 1342-016 
covering an area of 80" x 80". A large number of companion 
objects can be seen which, if at the redshift of the quasar, would 
appear consistent with a moderately rich cluster environment. 

of powerful AGN. In fact, at first sight, the il-band images 
presented in Appendix A and in McLure et al. (1999) show a 
relatively low occurrence of morphological disturbance, with 
only three objects (1012+008, 2349-014 & 0157+001) un- 
dergoing obvious, large-scale, interactions. However, as can 
be seen from the model-subtracted images, the removal of 
the best-fitting axisymmetric model for each host galaxy 
docs reveal the presence of peculiarities such as excess flux, 
tidal tails and close companions, at lower surface-brightness 
levels. It is tempting to conclude from this that many of 
these AGN may have been triggered into action by recent 
interaction with companion objects. However, the true sig- 
nificance of the levels of interaction seen in the AGN host 
galaxies can only be judged by comparison with an equally- 
detailed investigation of the morphologies of a control sam- 
ple of comparable, inactive galaxies. Here we attempt to 



quantify the significance of such peculiarities by constructing 
a control sample based on the twenty bright cluster galax- 
ies taken from the HST archive discussed in the previous 
subsection. 

With the range of redshifts covered by the five Abell 
clusters it was possible to construct 19 object pairs, each 
consisting of one of our AGN sample paired with a suit- 
able control, drawn from the 20 cluster galaxies. Model- 
subtracted images of each of the galaxy pairs were then 
assessed in a blind test for the occurrence of four possible 
indicators of disturbance or interaction, i.e: 

• Residual tidal/spiral arm features, 

• Residual asymmetric flux, 

• Residual symmetric fiux, 

• Companion object inside 10" radius, 

where neither the identity of the AGN, nor which object in 
each pair was active, was known to the decision maker. The 
results of this process arc shown in Table 7. 

The one clear result of this test is that, taken as a 
group, the active galaxies display a greater occurrence of 
tidal/spiral arm features in their model-subtracted images 
than the control sample. There is no significant difference 
in the other three indicators. It is noteworthy however, that 
of the six occurrences of tidal/spiral arms, only one of these 
(2349—014) is a radio-loud object. To test whether or not 
this indicates an inherent difference between the radio-loud 
and radio-quiet AGN, the blind test was repeated for all 33 
HST objects, without control galaxies, but still keeping the 
identity of each AGN secret, in order to prevent any subcon- 
scious bias. The results of this test are presented in Table 
8. 

Several features of the results presented in Table 8 merit 
comment. First, it can be seen that there is no significant 
difference between the RG and RLQ sub-samples, consistent 
with all the other results presented in this paper. Second, it 
would initially appear that there is a significant difference 
between the radio-loud objects and the RQQ sub-sample in 
terms of the occurrence of tidal/spiral arm features. How- 
ever, an examination of which of the RQQs are contributing 
to this result shows that four of them are the objects which 
have combined disc/bulge model fits. Therefore, if the com- 
parison is restricted to include only those RQQs with single- 
component elliptical host galaxies, this apparent difference 
disappears. 
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Table 7. The results of the blind test to determine any diflference in the prevalence of features indicative of interaction or disturbance 
in the AGN sample (active) compared to inactive galaxies, using 20 BCG as an inactive control sample (see text for details). 



Object Type Tidal/Spiral Asymmetric Flux Symmetric Flux Companion 
Active 5 8 5 7 



Inactive 10 7 9 



Table 8. The results of the blind test to determine if any differences exist in the occurrence of features indicative of interaction or 
disturbance between the AGN sub-samples (see text for details). 
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The results of these blind comparison tests thus support 
two separate conclusions. First, the finding that there are no 
detectable differences between the elliptical AGN host galax- 
ies and the BCG control sample provides further evidence 
that these two classes of object are directly comparable. Sec- 
ond, the apparently more frequent occurrence of tidal/spiral 
arm features in the RQQ sub-sample is only a manifestation 
of the previously-determined fact that a minority of RQQ 
host galaxies, two of which are not actually true quasars, 
have substantial disc components, and thus does not pro- 
vide evidence for an inherent difference in interaction rates 
between radio-quiet and radio-loud objects. 

Perhaps one should not be surprised by the lack of evi- 
dence for spectacular large-scale ULIRG-like interactions in 
our quasar sample, given that as described below, the acti- 
vation of a quasar requires the delivery of only ~ 1 M0yr~^ 
to the central black hole. 



8 THE AGN-HOST CONNECTION 

As discussed by McLure et al. (1999), our finding that the 
hosts of quasars are almost all spheroidal galaxies allows ex- 
ploitation of the now well-established correlation between 
black-hole mass and spheroid mass to estimate the masses 
of the central black holes in these objects, completely in- 
dependent of any of their observed nuclear properties. We 
can then use this information to explore which, if any, of 
the observed properties of the AGN are linked to black-hole 
mass. 

In this section we revisit this issue for a number of rea- 
sons. First, of course, our matched sub-samples are now com- 
plete. Second, in the intervening two years the application of 
3- integral models has resulted in a revising down of the con- 
stant of proportionality between Alth and Msph from 0.006 
(Magorrian et al. 1998) to ~ 0.003 - 0.001 (Kormendy & 
Gebhardt 2001). Third, independent estimates of the black- 
hole mass in a substantial number of our AGN have now 
been made on the basis of nuclear Hp line-width by McLure 
& Dunlop (2001). Fourth, the new VLA data presented in 
Section 4 mean that we can now explore the relationship 



between black-hole mass and radio luminosity without be- 
ing hampered by a large number of upper limits within the 
'radio-quiet' subsample. 

8.1 The Black-Hole Spheroid Connection 

As summarized by Kormendy & Gebhardt (2001), 'reliable' 
black-hole mass estimates are now available for at least 37 
nearby galaxies. As a result of this now substantial sample, 
clear correlations have been uncovered between black-hole 
mass and host-spheroid luminosity, and between black-hole 
mass and host stellar velocity dispersion. It has been claimed 
that the latter is the tighter correlation, indicating that the 
basic physical link is with spheroid mass. Since in the present 
study we do not possess measurements of host stellar veloc- 
ity dispersions, we obviously have to estimate host spheroid 
mass on the basis of the spheroid luminosity returned by our 
2-D modelling. However, a recent re-analysis of the bulge 
luminosities of low-redshift galaxies indicates that the re- 
lationship between bulge luminosity and black-hole mass is 
in fact just as tight as that between black-hole mass and 
central galaxy velocity dispersion (McLure & Dunlop 2002). 

To estimate black-hole masses on the basis of host 
spheroid luminosities we have therefore adopted the fol- 
lowing relations, while acknowledging the significant scatter 
present in both. 

Msph = 0.00123L^-^^ (7) 

Mbh = 0.0025M,pfe (8) 

consistent with the mass-to- light ratio relation for ellipticals 
determined by Jorgensen, Franx & Kjaergaaxd (1996) and 
the results of Gebhardt et al. (2000). We note that while 
these conversion factors arc certainly also consistent with 
those quoted in the recent review by Kormendy & Gebhardt 
(2001), there is now a growing body of evidence that the 
true conversion factor between Mgph and M^^ may be a 
factor ~ 2 lower than that adopted here (Merritt & Ferrarese 
2001, McLure & Dunlop 2002). If one wishes to adopt this 
lower factor (i.e. Mi,h = 0.0013M3p;j) then the estimates of 
black hole mass utilised in the analyses that follow in the 
remainder of this paper can simply be divided by a factor 
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of 2. Obviously none of the arguments concerning relative 
masses are affected by this uncertainty, and indeed none of 
the physical arguments are significantly altered by such a 
further modest reduction in the black-hole mass estimates. 

The results of applying these relations to calculate in- 
ferred bulge and black-hole masses are presented in Table 
9. Also given in Table 9 (column 4) are black-hole masses 
derived for the majority of the quasars in these samples by 
McLure & Dunlop (2001) on the basis of assuming that the 
iif/3-producing broad-linc clouds are in Keplerian orbits. 

The new spheroid-based black-hole estimates are com- 
pared with the J7/3-derived values in Fig 13. This diagram 
shows that, with few exceptions, there is very good agree- 
ment (~ 0.2 dcx) between the two completely-independent 
estimates of black-hole mass. Indeed, given the well-known 
uncertainties in both approaches, the level of agreement 
demonstrated by Fig 13 can be viewed as adding consid- 
erable credence to the estimated values of black-hole mass, 
and indeed also provide support for the basic premise of 
a gravitationally-bound broad-line region as assumed by 
McLure & Dunlop (2001, 2002). The fact that four objects 
(three visible on the diagram) have H^-hased black-hole 
masses which lie well below the corresponding host-spheroid 
derived values is not really surprising since the Hp line- width 
analysis is clearly capable of yielding a severe under-estimate 
of black-hole mass for any object whose broad-line region is 
orbiting close to the plane of the sky (but should be in- 
capable of yielding a comparably serious over-estimate of 
black-hole mass). Thus for 2247+140, 0736-^017, 0157-^001, 
and 1012-1-008 there is a good reason for trusting the black- 
hole mass estimates of ~ 10^ M© given in column 3 of Table 
9, more than the lower values yielded by the H/S analysis in 
column 4. 

Accepting this, and excluding the two low-luminosity 
RQQs, the third column of Table 9 shows that a black hole of 
minimum mass Mth > 5 x 10* Mq appears to be required to 
produce a rmclear luminosity corresponding to Mb < —23. 

Beyond confirming the basic plausibility of the black- 
hole mass estimates, the most interesting result contained in 
Fig 13 is that, as suggested by the initial results of McLure 
et al. (1999), the black- hole masses in the radio- loud quasars 
are systematically larger than in their radio-quiet counter- 
parts. In terms of summary statistics of central tendency, 
this difference, although clear, might not appear very dra- 
matic - specifically, excluding the two low-luminosity RQQs: 

{Mbh/l(f Mo) = 1.67 ± 0.36 median = 1.3 (KG) 

{Mbh/lO^ Mq) = 1.64 ± 0.22 median = 1.4 (RLQ) 
{Mbh/W'-' Mq) = 1.05 ± 0.24 median = 0.8 (RQQ) 

However, the difference becomes more striking when 
one notices that there is an apparent threshold mass of 
~ 10® Mq above which lie 8/10 RGs and 10/10 RLQs, but 
below which lie 9/13 of the RQQs. In other words, using 
the black-hole mass estimation procedure outlined above, 
at least within this particular sample it appears that a 
black-hole mass > lO'^ Mq is a necessary (albeit perhaps 
not sufficient) condition for the production of a powerful 
radio source. As measured via the KS test, this difference 
between the black-holo mass distributions of the radio-loud 
and radio-quiet subsamples is significant at the 3-a level, a 
significance which must be taken seriously given the almost 
perfect matching of the distributions of nuclear optical out- 
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Figure 13. A comparison between the black-hole masses pre- 
dicted from host-galaxy spheroid luminosity and those deter- 
mined from H/B line-width by McLure & Dunlop (2001). The 
shaded area illustrates that there is a region in which both ap- 
proaches agree that Mi,/^ > 10®'* Mq, and that this region con- 
tains all except 2 of the RLQs (open circles), but excludes all 
except 2 of the RQQs (filled circles). 



put demonstrated in Fig 4. Further support for its reality 
comes from the H(5 study by McLure & Dunlop (2001) al- 
ready referred to above, in which it was found that 11/13 
RLQs had Mm, > 10* ** Mq while only 4/18 RQQs had 
black-hole masses above this threshold. A similar conclu- 
sion was also reached by Laor (2000) from a study of the 
virial masses of PG quasars. 

An inevitable consequence of the matching of nuclear 
optical luminosities demonstrated in Fig 4, and the sys- 
tematic offset in black-hole mass described above is that, 
on average, the black holes at the heart of the RQQs in 
our sample must, on average, be emitting more efficiently 
than their radio-loud counterparts. This is quantified in the 
last two columns of Table 9 which give estimated Eddington 
luminosity (converted to i?-band absolute magnitude) and 
then emitting efficiency L/LEdd- The potential J?— band Ed- 
dington luminosity of each quasar was derived in the follow- 
ing manner. Initially, the host-galaxy bulge luminosities pro- 
vided by the two-dimensional modelling were converted into 
a spheroidal mass estimate using the J0rgensen et al. mass- 
to-light ratio (Eqn 7). The linear relation between black- 
hole mass and spheroidal mass (Eqn 8) was then used to 
estimate the central black-hole mass and resulting bolomet- 
ric Eddington luminosity of each quasar. To calculate the 
corresponding i?— band magnitudes it was assumed that, 
on average, the bolometric luminosity of a quasar can be 
estimated via Ltoi — IOAL5100 (Laor ct al. 1997), where 
L5100 is the absolute monochromatic luminosity at 5100 A. 
The resulting estimate of the absolute luminosity at 5100 
A was then transformed to the central wavelength of the 
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Table 9. The results of estimating black-hole mass from the i?-band luminosity of the host spheroid. Columns two and three list the 
predicted galaxy spheroid mass and central black-hole mass respectively. For comparison with the results given in Column 3, Column 
4 gives the black-hole mass estimates as derived from the Hp line width by McLurc & Dunlop (2001). Colurrm 5 lists the predicted 
absolute i?-band Eddington luminosity of the black-hole, calculated according to the prescription given in Section 8.1. Column 6 gives 
the ratio of the predicted Eddington luminosity to the best-fitting nuclear model component. 



Source 


rusph/ J^u ivi0 


nibh/ J^U iVl0 


Hlbh— H/3/ J^U iVl0 




T IT,, 


RG 












0230-027 


4.8 


1.2 








0307+169 


7.9 


2.0 








0345+337 


2.6 


0.7 








0917+459 


10.6 


2.6 








0958+291 


3.8 


0.9 








1215+013 


2.1 


0.5 








1215-033 


3.8 


0.9 








1330+022 


5.8 


1.4 








1342-016 


16.7 


4.2 








91 A^ -1-97Q 


y.o 


9 'i 
z.o 








RLQ 












0137+012 


8.7 


2.2 


1.1 


-28.2 


0.013 


0736+017 


5.0 


1.3 


0.3 


-27.6 


0.035 


1004+130 


9.4 


2.3 


3.7 


-28.3 


0.089 


1020-103 


3.8 


1.0 


0.7 


-27.4 


0.028 


1217+023 


5.9 


1.5 


0.8 


-27.8 


0.043 


2135-147 


3.9 


1.0 


2.6 


-27.4 


0.049 


2141+175 


4.2 


1.1 


1.7 


-27.5 


0.061 


2247+140 


6.5 


1.6 


0.1 


-27.9 


0.022 


2349-014 


12.2 


3.1 


1.8 


-28.6 


0.014 


2355-082 


5.3 


1.3 


0.7 


-27.7 


0.014 


RQQ 












0052+251 


2.3 


0.6 


0.6 


-26.8 


0.107 


0054+144 


5.2 


1.3 


2.5 


-27.7 


0.052 


0157+001 


12.2 


3.1 


0.2 


-28.6 


0.011 


0204+292 


3.6 


0.9 


0.6 


-27.3 


0.020 


0244+194 


1.9 


0.5 


0.3 


-26.6 


0.046 


0257+024 


0.4 


0.1 




-24.9 


0.008 


0923+201 


3.4 


0.8 


2.6 


-27.2 


0.088 


0953+414 


2.1 


0.5 


0.7 


-26.6 


0.355 


1012+008 


6.4 


1.6 


0.02 


-27.9 


0.026 


1549+203 


1.0 


0.3 




-25.9 


0.169 


1635+119 


3.2 


0.7 


0.4 


-26.9 


0.007 


2215-037 


4.7 


1.2 




-27.6 


0.010 


2344+184 


1.3 


0.3 




-26.8 


0.002 



R— filter (6500 A) assuming a power-law quasar spectrum 
(a — 0.2). Using tliis prescription, a quasar with a central 
black-hole mass of 10® Mq will have an Eddington-limit 
absolute il- band magnitude of Mu.[Edd) = —27.5. 



Observed nuclear Mb is plotted against potential 
Eddington-limited Mn in Fig 14. As can be seen from this 
figure, the maximum luminosity produced by any quasar 
is within a factor of two of the predicted Eddington limit, 
while the majority appear to be radiating at ~ 1 — 10% of 
the Eddington luminosity. Taken as a group there is a sug- 
gestion that the RQQs are radiating at a higher percentage 
of the Eddington limit, but the difference is not formally 
significant (KS, p = 0.29). 



8.2 The Black-Hole— Radio Connection 

When the matched samples of RLQs and RQQs were first 
defined, all that was known about the radio properties of 
the RQQs was that their 5 GHz luminosities were lower 
than PsGHz ^ 10^^ WHz^^sr"^ (Dunlop ct al. 1993). How- 
ever, with the completion of the VLA detection programme 
described in Section 4 we can now explore how the radio 
properties of all the AGN in our combined sample relate to 
their estimated black-hole masses. 

In Fig 15 total radio luminosity, P^ghz, is plotted 
against estimated black-hole mass for all the objects in our 
combined sample (with the single exception of 1549+203, 
which was not observed with the VLA to the same depth 
as the other objects). Also included in this diagram are the 
appropriate points for 8 nearby galaxies taJien from Prances- 
chini et al. (1998), with the Milky Way and M87 highUghted. 
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Several significant trends are evident in this diagram. 
First it can be seen that, within the RQQ sample, radio lu- 
minosity correlates, albeit weakly, with black-hole mass, al- 
though our sample doesn't span a wide dynamic range and 
the correlation is not formally significant. Second, and per- 
haps more interestingly, the RQQs are consistent with the 
locus described by the nearby galaxies (for which, of course, 
the black-hole masses are based directly on stellar dynam- 
ics). This simultaneously reinforces the plausibility of our 
black-hole mass estimates, and suggests a common physical 
origin for the radio emission from RQQs and nearby opti- 
cally more-quiescent bulges. Whatever this physical mech- 
anism is, it must explain the observed strong dependence 
of minimum radio output on black-hole mass, in which the 
power- law index 7 is at least 2 (adopting Pradio oc M^f^. As 
noted by Franceschini et al. (1998), 7 = 2.2 is the prediction 
of ADAF models, while 7 ~ 2 is a generic prediction of any 
mechanism in which radio output is proportional to the area 
of the accretion disc. 

Our results therefore provide further support for pre- 
vious claims that there is a minimum radio output from 
black holes and that this is a strong function of mass. How- 
ever, while this obviously implies that the most massive 
black holes can never be very radio-quiet. Fig 15 graphically 
demonstrates that the genuine radio-loud objects lie on a 
quite distinct relation. Interestingly, our radio-loud sample 
also displays an (albeit weak) internal correlation between 
PsGHz and estimated black-hole mass. Either this relation- 
ship must be much steeper (P oc M^^) than the minimum 
radio-output relation, or, as illustrated in Fig 15, the black- 
hole mass dependence of the upper envelope may be con- 
sistent with that of the lower envelope but offset by several 
orders of magnitude. 

In fact the findings of Lacy et al. (2001) strongly favour 
the latter option. They report a relationship between the 
logarithm of radio power and black-hole mass which has a 
slope of 1.4 ± 0.2. However, inspection of their Fig 2 shows 
that their data are also consistent with an upper and lower 
envelope with slope ~ 2.5 as illustrated here in Fig 15. 

We note here that it has recently been claimed by Osh- 
lack, Webster & Whiting (2001) that the radio-loud Seyfert 
gala^xy PKS 2004—447 contradicts the existence of such an 
envelope, and shows that relatively low-mass black holes 
can produce a powerful radio source. However, this claim 
is based on their black-hole mass estimate for this object of 
~ 5 X 10^ M0 and it seems almost certain that this value, de- 
rived from the relatively narrow emission lines displayed by 
PKS 2004—447, is a serious under-estimate of the true mass. 
The reason for this is that since this object displays optically 
variability and has a compact radio source, it seems very 
likely that its orbiting broad-line region lies close to plane 
of the sky. In this situation (as found in our own sample 
for, e.g. 1012-1-008 - see Table 9) the black-hole mass de- 
rived from H/j emission-line width can under-estimate the 
true value by a factor > 100. In Fig 15 we have plotted the 
position of PKS 2004—447 after re-calculating its black-hole 
mass on the assumption that the orbital plane of its broad- 
line region is oriented within ~ 5 degrees of the plane of 
the sky (sec McLurc & Dunlop 2001). This increases the es- 
timated black-hole mass to ~ 4 x 10*^ M©, and moves PKS 
2004—447 into complete consistency with the upper enve- 
lope delineated by the other available data. We have not 
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Figure 14. The observed absolute magnitude of the nuclear 
component in each quasar plotted against the absolute magnitude 
which is predicted by assuming that each quasar contains a black 
hole of mass mi,h = 0.0025mspheroid, and that the black hole 
is emitting at the Eddington luminosity (RLQs = open circles, 
RQQs = filled circles). The solid line shows where the quasars 
should lie if they were all radiating at their respective Edding- 
ton luminosities, while the dashed line indicates 10% of predicted 
Eddington luminosity, and the dot-dashed line indicates 1% of 
predicted Eddington luminosity. The nuclear components of the 
radio galaxies are not plotted because all the evidence suggests 
they are substantially obscured by dust - see Section 6.7. 

corrected the observed radio emission for dopplcr boosting - 
if we did then the this object would simply move downwards 
in Fig 15. 

Thus both our own results and other currently available 
data support the existence of both a lower and an upper 
envelope to the radio luminosity which can be produced by 
a black hole of given mass. Indeed Fig 15 serves to provide 
an elegant explanation not only of our own findings, but also 
of several other well-known results in the literature. 

First, because we selected our radio-loud objects on 
the basis that their radio luminosities were greater than 
/opio(P/WHz~^sr~^) = 24, the form of the upper envelope 
in these figures can be seen to be perfectly consistent with 
our finding that virtually all our radio-loud objects have 
Mbh > W'' Mq. In the context of this diagram this finding 
can be seen to be a result of our radio-selection criterion 
rather than the existence of some critical black-hole mass 
for the production of relativistic jets within a galaxy mass 
halo. 

Second, the steep dependence of this apparent upper 
envelope on black hole mass means that radio sources up to 
/opio(P/WHz~^sr~^) = 27 — 28, consistent with the most 
luminous SCR galaxies even at z = 1 can be produced by 
black holes (and hence hosts) only a factor of 2-3 more mas- 
sive, and most likely will be due to the steep high-mass slope 
of the black-hole mass function. This provides a natural ex- 
planation of why radio galaxies spanning the top two orders 
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of magnitude in radio power all appear to be giant ellipti- 
cals with stellar masses which agree to within a factor of ~ 2 
(Jarvis et al. 2001). 

Third, from Fig 15 it can be seen that the properties 
of AGN hosts will inevitably become more mixed if the 
radio- luminosity threshold is move down to, say, Psghz — 
10^^ WHz^^sr"^. Given the lower boundary in Fig 15, ra- 
dio sources in this power range can still be produced by 
black holes as massive as 10'^ Mq living in giant ellipticals, 
but it is also clear that black holes with masses as low as 
IQ^-s are capable of producing detectable radio sources 
at this level. As demonstrated by the properties of the lower- 
luminosity objects in our own sample, and by the studies of 
'normal' galaxies at low redshift, such black holes can be 
housed in disc-dominated galaxies as well as discless lower- 
mass spheroids. 



pks 2004-447 ^- / 




8.3 The origin of radio loudness 

While the broad dependence of radio luminosity on black- 
hole mass may explain the relative numbers of radio-quiet 
and radio-loud objects, it is also clear from Fig 15 that black 
hole mass cannot be invoked to explain the range of radio 
luminosities displayed by quasars of very similar optical lu- 
minosity. Specifically it is clear from this figure that objects 
powered by equally massive black holes can differ in radio 
luminosity by ~ 4 orders of magnitude. 

Interestingly Fig 15 also allows us to rule out the pos- 
sibility that the radio output of a black hole of given mass 
is driven simply by fuel supply. This is because, while it has 
long been known that RQQs are not radio silent, it can be 
seen from Fig 15 that several of the RQQs in our sample arc 
in fact as radio-silent as they could possibly be, given their 
estimated black-hole masses. In other words at least some 
of the black holes in RQQs produce as little radio output as 
their counterparts in completely quiescent nearby galaxies, 
despite the fact that black holes in quasars are clearly in re- 
ceipt of sufficient fuel to produce powerful optical emission. 

Thus, while the results of this study cannot by them- 
selves provide a definitive answer to the long standing ques- 
tion of the origin of radio loudness, they can be used to focus 
the argument by excluding several possible explanations. In 
summary, our results can be used to exclude host galaxy 
morphology, black-hole mass, or black hole fueling rate as 
the primary physical cause of radio loudness. 

We are thus forced to the conclusion that the produc- 
tion of powerful relativistic radio jets is driven by some other 
property of the black hole itself, the most likely candidate 
being black-hole spin (e.g. Blaridford 2000, Wilson & Col- 
bert 1995). Indeed, following the suggestion by Blandford 
(2000), we speculate that the evolutionary track followed by 
an activated black hole may manifest itself as a near verti- 
cal descent in Fig 15, with an active rapidly spinning hole 
appearing first as a powerful radio source close to the upper 
envelope in this diagram, and then descending towards the 
lower envelope as the spin energy of the hole is exhausted. 



LogioCmbh/M©) 

Figure 15. Total radio luminosity P^q'^^ versus black-hole mass 
showing the data on low-redshlft 'normal' galaxies from Prances- 
chini et al. (1998) (as squares, with the Milky Way and M87 high- 
lighted), and the AGN from the HST sample (RGs = crosses, 
RLQs = open circles, RQQs = filled circles). The lines which 
bracket the data are relations of the form P^CHz oc rn^'^ offset 
by 2.5 orders of magnitude from each other. The solid line there- 
fore indicates the apparent existence of a lower limit to the radio 
output of a black hole which scales as rnf^'^ while the dot-dash line 
indicates the apparent upper limit required by our data. While our 
own data are insufficient to constrain the slope of such an upper 
envelope, adoption of the same slope as the lower envelope (as 
done here) is consistent with the data of Lacy et al. (2001) and is 
further reinforced by the results of McLure & Dunlop (in prepara- 
tion). PKS 2004—447 has been included in this plot because it has 
recently been claimed that this object contradicts the existence 
of any such mass-dependent upper envelope (Oshlack, Webster & 
Whiting 2001). The position of this compact and variable source 
on this diagram has been here re-calculated on the assumption 
that its radio axis lies with ~ 5 degrees of the line of sight, but 
without maJjing any corrections for relativistic beaming of its ra- 
dio output (see text for details). As can be seen, the resulting 
increase in estimated black hole mass from ~ 5 x 10® Mq to 
mih ~ 4 X 10* M0 is sufficient to bring it into consistency with 
the upper envelope suggested by the other data. 

9 COSMOLOGICAL IMPLICATIONS 

9.1 Relative numbers of radio- loud and 
radio-quiet quasars 

One interesting aspect of the different black-hole mass 
thresholds for RQQs and RLQs uncovered by this study 
is that this difference of a factor of 2 in minimum black- 
hole mass provides a natural explanation of why radio-quiet 
quasars outnumber their radio-loud counterparts by a fac- 
tor of ~ 10 (and why this factor may reduce with increasing 
optical luminosity; Hooper et al. 1995, Goldschmidt et al. 
1999). 

This is because if our black-hole mass thresholds of 
Mbh > 5 X 10* Mq for RQQs, and Mth > 1 x 10^ for 
radio-loud objects are converted back into absolute spheroid 
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magnitudes in the if -band (using the absolute _R-band 
magnitudes given in Table 4, and then converting into K 
adopting R — K = 2.5 as justified by the results in Sec- 
tion 6.6) we find that RQQs can be hosted by spheroids 
with Mk < —25.3 while radio-loud objects require a host 
spheroid with Mk < —26. With reference to the Jf-band 
luminosity function these absolute magnitude limits corre- 
spond to L > 1.5L* and L > 3L* respectively (Gardner et 
al. 1997, Szolkoly et al. 1998, Kochanek et al. 2001) which 
loads to the prediction that, in the present day universe there 
are 1 x 10"'* galaxies per Mpc^ capable of hosting an RQQ 
brighter than My = —23, but only 1 x 10~® galaxies per 
Mpc^ arc capable of hosting an RLQ. This means that, as- 
suming low-redshift quasars arise from a raudomly-triggered 
subset of massive black-holes in spheroids, the factor of two 
difference in black-hole mass threshold translates into a fac- 
tor of 10 difference in the number density of RQQs and RLQs 
due simply to the steepness of the bright end of the spheroid 
luminosity function. 

9.2 Black-hole activation fraction at z ~ 

We can go further and use the above numbers to estimate 
the activation fraction for massive black holes in the local 
universe. First we note that the radio-loud objects in our 
sample have P2.7GHZ > 1 x lO'^^ WHz^^sr^^ (assuming a 
spectral index a ~ 0.8). Reference to the local radio lumi- 
nosity function in Dunlop & Peacock (1990), or to the more 
recent determination at 1.4 GHz from the 2dF by Sadler 
et al. (2001), indicates that the present-day number density 
of radio sources above this threshold is ~ 1 x 10~* Mpc~^. 
Thus, in the present-day universe we find that 1 in every 
1000 black holes more massive than 1 x lO'^ M0 (or equiva- 
lently spheroids more luminous than 31/*) is radio active in 
the present-day universe. 

Turning to the RQQ population, the present-day rmm- 
ber density of quasars with My < — 23 can be estimated by 
extrapolating the QSO OLF deduced at « ~ 0.4 by Boyle 
et al. (2000) back to 2 = assuming luminosity evolution 
oc (1 -I- z)^ . This produces an estimated number density of 
~ 5 x 10^^ Mpc^"', which is at least consistent with the di- 
rect determination of the low- 2; OLF attempted by Londish, 
Boyle & Schade (2000). Boosting this number by a factor 
of two to allow for obscured quasars (i.e. adopting an open- 
ing angle of ~ 45° for consistency with the correction factor 
which appears to apply to the radio-loud population) then 
leads to the conclusion that 1 in every 1000 black holes more 
massive than 5 x 10* M© (or equivalently spheroids more 
luminous than 1.5L*) is producing quasar-level optical ac- 
tivity. 

Thus both these comparisons yield the same result, 
namely that in the present-day {z < 0.1) universe 1 in ev- 
ery 1000 massive black holes is active, and that the reason 
RQQs outnumber RLQs by a factor of 10 is a direct result 
of the fact that their respective minimum black-hole masses 
differ by a factor of two. 

We note here that Wisotzki, Kuhlbrodt & Jahnke 
(2001) have recently concluded that the present-day black- 
hole activation fraction is 1 in 10000, or 1 in 5000 if the 
quasar LF is boosted by a factor of two as assumed above. 
However, reference to their Fig 2 shows that this factor re- 
sults from adopting the spheroid luminosity function of Lin 



et al. (1999). If instead Wisotzki et al. were to adopt the 
Kochanek et al. (2001) luminosity function then it is clear 
that (at least for the mass range probed by the present 
study) Wisotzki et al.'s analysis would also yield an acti- 
vation fraction of 0.1%. 



9.3 Black-hole activation fraction at z ~ 2 — 3 

At the peak epoch of quasar activity between z ~ 2 and 
2 ~ 3, the co-moving number density of powerful radio 
sources with P2.7GHZ > 1 x 10^'' WHz"*sr^* has risen to 
~ 1 X 10"'' Mpc"^ (Dunlop & Peacock 1990) while that 
of optically-selected quasars brighter than My = —23 is 
~ 5 X 10"® Mpc"^ (Warren, Hewett & Osmer 1995, Boyle 
et al. 2000). Doubling the latter figure to make the same cor- 
rection for obscured RQQs as applied in the previous sub- 
section at low redshift, leads to the conclusion that, as at 
2 ~ 0, RQQs at 2 ~ 2.5 outnumber their radio-loud counter- 
parts by a factor of ~ 10. However, at z ~ 2.5 the co-moving 
number density of both classes of object is enhanced by a 
factor of 100, giving an activation fraction of massive black 
holes of ~ 10%, assuming that the entire massive black-hole 
population is in place by that epoch. 

9.4 The cosmological evolution of AGN 

A renewed attempt to model the cosmological evolution of 
AGN in the light of our host-galaxy results is obviously 
beyond the scope of this paper. However, if one assumes 
that the same mass thresholds uncovered in our low-redshift 
study apply at high redshift, then one simple interpretation 
of the dramatic (and very similar) evolution of radio-loud 
and radio-quiet AGN between « ~ 2.5 and the present-day 
is that it is primarily due to density evolution, with the ac- 
tive fraction of massive black holes dropping from 10% to 
0.1% over this period. There must of course also be an ele- 
ment of luminosity evolution because pure density evolution 
was excluded long ago as an acceptable representation of 
the evolving RLF or OLF. However, as suggested by Miller, 
Percival & Lambert (in preparation) if individual quasars 
produce declining light curves it is inevitable that a larger 
fraction will be observed closer to peak luminosity as the 
activation rate rises, and Fig 14 indicates that at least some 
of the host galaxies uncovered in the present study should 
be capable of producing quasars with My ~ —28 if observed 
while accreting at the Eddington limit. 

One prediction of this scenario is that, as a larger frac- 
tion of massive black holes become more efficiently fueled, 
the apparent mass difference between RLQ and RQQ hosts 
should grow with increasing redshift. The reason for this 
is that, if fueling rates are, on average, somewhat higher 
at high redshift we can expect optically selected QSOs of 
a given luminosity to be produced by black holes with, on 
average, lower masses than at low redshift. However, if the 
rribh > 10^ Mq mass threshold for powerful radio activity 
found in this study continues to apply at high redshift, then 
radio selection will continue to yield only the most mas- 
sive black holes residing (presumably) in the most massive 
spheroids. 

In fact, studies of the hosts of RLQs and RQQs out to 
2 ~ 2 are already uncovering evidence of just such a trend. 
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Specifically, Kukula et al. (2001) have used NICMOS on the 
HST to measure the rest-frame _R-band luminosities of the 
hosts of matched samples of RLQs and RQQs at z ~ 1 and 
a ~ 2. Kukula et al. find that the hosts of RLQs are essen- 
tially unchanged in mass over this rodshift range, suggesting 
that the minimum mass threshold for powerful radio emis- 
sion indicated by the upper locus in Fig 15 applies at all 
redshifts and is therefore of physical significance. However, 
they also find that the average ratio of RLQ:RQQ host lu- 
minosity rises from the value of 1.5 found here at z ~ 0.2, 
to ~ 2 at 2; ~ 1, and ~ 3 at 2; ~ 2, suggesting a progres- 
sive drop in the average mass of RQQ hosts with increasing 
redshift, a finding supported by the results of other recent 
studies (Ridgway et al. 2000, Rix et al. 2000). 



10 CONCLUSIONS 

In this paper we have reported the extensive results which 
follow from completion of our HST imaging programme of 
the host galaxies of low-redshift radio-quiet quasars, radio- 
loud quasars, and radio galaxies. This paper represents the 
completion of a programme which has spanned most of the 
last decade, commencing with the deep infrared imaging of 
the same sample of objects by Dunlop et al. (1993). The 
depth and quality of our HST data have now allowed us to 
determine accurately all the basic structural parameters of 
the hosts of these 3 classes of powerful AGN. Because of this, 
and because of the wealth of data now available for this sam- 
ple at other wavelengths (including the new deep VLA ob- 
servations also reported here), we have been able to address 
a wide range of issues of which are hopefully of interest to 
workers in several different areas of extra-galactic research. 
Consequently we conclude with a detailed summary of the 
main conclusions of this study, structured to assist the in- 
terested reader in moving directly to the sections which may 
be of most relevance to their own work. 

10.1 Results from analysis of the HST images 

Prom the detailed 2-dimensional modelling of the new HST 
images presented here in combination with the data reported 
by McLure et al. (1999) we find that (as detailed in Sections 
5 and 6): 

• All except the two least luminous RQQs in the sample 
have bulge-dominated hosts, and in only 2 of the remaining 
31 objects (the RQQs 0052-^251 and 0157-1-001) can we find 
any evidence for a significant disc component. Thus both 
radio-loud and radio-quiet quasars live, almost universally, 
in elliptical galaxies. 

• The hosts of all throe classes of powerful AGN are lu- 
minous galaxies with L > L* , and almost always L > 2L* . 
The average luminosities of the RG and RLQ hosts are es- 
sentially identical (consistent with radio- loud unification), 
and equivalent to 41/*. The average luminosity of the RQQ 
hosts is somewhat smaller, equivalent to ~ 3L* . 

• The hosts of all three classes of powerful AGN are large 
galaxies with typical half-light radii ri/2 — 10 kpc. As with 
the luminosities, the average values of host-galaxy scale- 
length for the radio-loud sub-samples are indistinguishable, 
while the hosts of the RQQs are typically smaller by a mod- 
est factor (~ 1.5). 



• The hosts of all three classes of powerful AGN display 
a surface-brightness scalelength (Kormendy) relation iden- 
tical (in both slope and normalization) to that displayed by 
inactive massive ellipticals predominantly found in clusters. 

• The hosts of all three classes of powerful AGN display 
a distribution of axial ratio which is indistinguishable to 
that which has been long established for normal elliptical 
galaxies. 

10.2 Results incorporating infrared images 

From a joint analysis of the HST optical images with existing 
K-hand images of the same sample (Dunlop et al. 1993, 
Taylor et al. 1996, McLure, Dunlop & Kukula 2000) we find 
that (as detailed in subsections 6.6 and 6.7): 

• The hosts of all three classes of powerful AGN are 

red galaxies, with typical rest-frame optical-infrared colours 
R — K ~ 2.5. Consistent with the results of Nolan et al. 
(2001), these colours are as expected from an evolved stel- 
lar population of age 10 — 13 Gyr, indicating that, as for 
quiescent massive ellipticals, the stellar mass of an AGN 
host galaxy is dominated by a well-evolved stellar popula- 
tion formed at high redshift (« > 2). 

• The optical-infrared colours of the RQQ and RLQ nu- 
clei are statistically indistinguishable, and lie in the range 
R — K ~ 2 — 4. For the RG nuclei, the relation between 
R — K colour and R-haiiA luminosity is as expected under 
the assumption that they all contain obscured nuclei with 
intrinsic optical luminosities comparable to the quasars, red- 
dened by dust with kx cc A'^ ''^. This result strongly favours 
unification of RGs and RLQs via orientation rather than 
time. 

10.3 Relation of quasar hosts to normal galaxies 

Prom a comparison of our results with those derived from 
recent studies of the 'normal' massive galaxy population, 
and of low-rcdshift ULIRGs and Seyferts, we find that (as 
detailed in Section 7): 

• The bulge:disc ratio of AGN hosts is a function of AGN 
luminosity, with disc-dominated hosts dying out above the 
traditional quasar:Seyfert boundary at Mv{nuc) ~ —23. 

• In contrast to ULIRGs, quasar hosts lie in a region 
of the fundamental plane (as judged from the photomet- 
ric projection of the FP offered by the Kormendy relation) 
which has been shown to be occupied by the most massive 
and apparently old population of ellipticals which display 
'boxy' isophotes and distinct kinematic cores. This provides 
a strong argument against the possibility that the r^'^'*-law 
luminosity profiles displayed by the vast majority of the 
quasar hosts in our sample are the result of recent merg- 
ers between massive disc galaxies. This result also offers lit- 
tle support for a strong evolutionary connection between 
ULIRGS and quasars in the low-rcdshift universe. 

• The basic structural properties of AGN hosts axe indis- 
tinguishable from those of inactive brightest cluster galaxies. 
Consistent with the environmental study of McLure & Dun- 
lop (2000), the largest host galaxies in our sample are of 
comparable mass to the BCGs found at the heart of Abell 
Class 1 or Class 2 clusters. However, the lack of any host 
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galaxies as massive as an Abell Class 4 cluster BCG is as 
expected given the relative rarity of such rich environments. 
All the evidence considered is consistent with the AGN hosts 
being drawn, essentially at random, from the present-day 
massive-elliptical galaxy population. 

• There is no statistically-significant evidence that AGN 
hosts display more signs of large (kpc) scale disturbance, or 
multiple nuclei than normal comparably-massive quiescent 
ellipticals. 

10.4 The black-hole spheroid connection 

Combining our HST data with our new VLA data, and util- 
ising recent results on the black-hole:spheroid mass relation 
in both quiescent and active galaxies, we find that (as de- 
tailed in Section 8): 

• Estimates of central black-hole mass for the quasars in 
our sample based on host-galaxy spheroid luminosity agree 
well with those derived from Hp emission-line width by 
McLure & Dunlop (2001). For three objects the latter tech- 
nique yields a much lower value, but this is as expected for 
a small subset of objects if the broad-line region has a pre- 
dominantly disc-like structure. 

• Based on the relation = 0.0025Msp;i , we find that 
all the quasars in our sample are powered by a black-hole of 
mass rribh > 5 x 10* Mq , but that the radio-loud objects lie 
above an even higher mass threshold, with mth > 10^ M©. 

• The most efficient RQQ in our sample appears to be 
radiating at ~ 35% of the Eddington limit, but the vast 
majority of these low-redshift AGN axe emitting at 1 — 10% 
of the Eddington limit. Based on these calculations, the most 
massive objects in this low-rcdshift sample appear capable 
(if fuelled at maximum efficiency) of producing quasars with 
Mv ^ —28, comparable to the most luminous objects known 
at high redshift. 

• The black-hole mass difference between RLQs and 
RQQs reflects the existence an apparent upper and lower 
limit to the radio output that can be produced by a black 
hole of a given mass. Both the upper and lower thresholds 
on radio luminosity appear to be a strong function of black 
mass (oc 'ml'/^^"')- At least some RQQs in our sample appear 
to be as radio- weak as 'normal' inactive comparably-massive 
spheroids found in the low-redshift universe, despite the fact 
that must be receiving sufficient fuel to power the optical 
quasar nucleus. 

• While there is a broad and clear trend for increasing 
radio luminosity with increasing black-hole mass, it is now 
cleax that host morphology, black-hole mass, and fuel supply 
can each be excluded as the primary physical explanation of 
why a subset of quasars arc up to 4 orders of magnitude 
more radio luminous than their 'radio-quiet' counterparts. 
Consequently we argue that black-hole spin is the most plau- 
sible (perhaps only) remaining feasible explanation for the 
production of a powerful radio source. 

10.5 Cosmological implications 

Considering the implications of our results for our under- 
standing of the nature and evolution of AGN populations as 
a function of redshift, we find that (as detailed in Section 
9): 



• The relative numbers of radio-loud and radio-quiet 
quasars can be naturally explained by the above-mentioned 
mass thresholds, combined with the form of the bright end 
of the elliptical galaxy luminosity function. 

• The activation fraction of massive black holes is ~ 0.1% 
in the present-day universe, rising to ~ 10% at « ~ 2 — 3, 
corresponding to the peak epoch of quasar activity in the 
universe. 

• The black-hole mass threshold for powerful radio activ- 
ity can explain the trend for a growing gap with increas- 
ing redshift between the masses of RLQ hosts and RQQ 
hosts uncovered by recent HST-based studies of high red- 
shift quasars. 
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APPENDIX A: NEW HOST-GALAXY IMAGES 
AND MODEL FITS 

The images, two-dimensional model fits, and model sub- 
tracted images for the 14 AGN for which the new HST data 
are reported in this paper are presented in this appendix 
in Figs Al to A14. A grcy-scalo/contour imago of the fi- 
nal reduced F675W -R-band image of each AGN is shown 
in the top-left panel (panel A) of each figure, covering a 
region of 12.5 x 12.5 arcsec centred on the target source. 
The surface-brightness level of the lowest contour is indi- 
cated in the top-right corner of this panel, and the grey-scale 
has been chosen to highlight structure close to this limit. 
Higher surface-brightness contours are spaced at intervals 
of 0.5 mag. arcsec"^, and have been super- imposed to em- 
phasize brighter structure in the centre of the galaxy/quasar. 
Panel B in each figure shows the best-fitting two-dimensional 
model, complete with the unresolved nuclear component (af- 
ter convolution with the empirical PSF described above), 
contoured in a manner identical to panel A. Panel C shows 
the best-fitting host galaxy as it would appear if the nuclear 
component were absent, while panel D is the residual image 
which results from subtraction of the full two-dimensional 
model (in panel B) from the raw i?-band image (in panel 
A), in order to highlight the presence of any morphological 
peculiarities such as tidal tails, interacting companion galax- 
ies, or secondary nuclei. Within each figure, all four panels 
are displayed using the same greyscale. 

Al Notes on individual objects 

Here we provide a brief discussion of each of the 14 new 

HST images, with reference to other recent HST and ground- 
based data. Comparably-detailed descriptions of the other 
19 objects in the sample can be found in McLure et al. (1999) 
and are not repeated here. Additional details on each object, 
along with a description of the main features of our K-band 
images, can be found in Dunlop et al. (1993) and Taylor et 
al. (1996). The results of off- nuclear optical spectroscopy and 
resulting spectral model-fitting for several of the objects in 
this sample can be found in Hughes et al. (2000) and Nolan 
et al. (2001). 

A 1.1 The Radio Galaxies 

0230-027 (PKS 0230-027, CD -050) 

The now _R-band HST imago shows this galaxy to be uniform 
and round with no sign of obvious distortion. The detail of 
the WF2 image shown in Fig Al reveals ~ 10 apparent com- 
panion objects lying in a roughly circular formation around 
the galaxy, at a radius of ~ 30 — * 60 kpc. The brightest of 
these companions lies at a projected distance of ~ 60 kpc to 
the NW. The full WF2 chip reveals there to be ~ 20 faint 
companion objects within a radius of ~ 200kpc. 

0307+169 (3C 079, 4C -hl6.07) 

The R-hsnd image of this source shown in Fig A2 reveals it 
to be a classic brightest-cluster galaxy (BCG). Two bright 
companion objects can be seen to the North and South with 
three or four accompanying tidal arm features emanating 
from the central galaxy. The suggestion that this source is 
in the process of undergoing merger activity is strengthened 



by the overlying contours which show three distinct cores. 
The image of the full WF2 chip reveals several more bright 
companions within a radius of < 100 kpc. 

This source has recently been imaged by the HST PC 
during the 3CR snapshot survey (Dc Koff ct al. 1996) 
through the F702W (wide R) filter. The 280-second inte- 
gration presented by De Koff et al. confirms the complex 
multiple structure of this object with the authors noting 
that the optical and radio axes are aligned to within 15 deg. 
McCarthy ct al. (1995) imaged this source both in the R- 
band and through an Ha emission-line filter, detecting a 
curving filament of extended Ha emission stretching some 
12" to the NW. 

1215-033 (PKS 1215-033) 

The -R-band image of this source shows it to be a uniform 
round galaxy with no obvious signs of interaction or distur- 
bance. There are two companion objects lying just on the 
edge of the detail shown in Fig A3, to the SW and ESE at 
projected distances of ~ 60 kpc. 

1215+013 (PKS 1215+013) 

The _R-band image of this source shown in Fig A4 reveals nu- 
merous faint companion objects around the central galaxy, 
the brightest lying ~ 30 kpc to the South. The overlying 
contours for this source suggest some sort of disturbance 
immediately to the North of the galaxy core. 

1330+022 (3C 287.1, 4C +02.36) 

Numerous companion objects are seen in the new _R-band 
image of this source shown in Fig A5, the brightest of which 
can be seen some 25 kpc to the NE. There is a linear tidal 
feature ~ 35 kpc to the SW. The overlying contours for this 
source show there to be an apparent second nucleus at less 
than 1" separation to the WNW 

This object was in the sample imaged in the i?-band 
during the 3CR snapshot survey (De Koff et al. 1996). The 
280-second image presented by De Koff et al. confirms the 
presence of the second rmclcus. 3C 287.1 has been shown 
to have a power-law X-ray spectrum by Crawford & Fabian 
(1995). 

1342-016 (PKS 1342-016, MRC 1342-016) 
The new HST image shows this galaxy to be large, luminous 
and uniform. A bright foreground star is present ~ 13" to 
the SWW, on the edge of the frame shown in Fig A6. No 
obvious companion objects are present in the detail shown 
in Fig A6 although the full WF2 chip image shown in Fig 
13 reveals there to be a large number of companion objects, 
consistent with a moderately rich cluster. 

A 1.2 The Radio-Loud Quasars 

1020-103 (PKS 1020-103, UT 1020-103) 
The i?-band image of this source shown in Fig A7 reveals 
it to have a comparatively small and faint host, with an 
obvious PA of ~ 110 degrees. There is a triangle of three 
faint companions directly to the South. This object was 
identified as an X-ray source in the ROSAT All-Sky Sur- 
vey (Brinkmann et al. 1997). 

1217+023 (PKS 1217+02, ON 029) 

The J?-band image of this relatively nuclear-dominated ob- 
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ject shown in Fig A8 clearly reveals an elliptical-looking host 
galaxy with an apparent PA of ~ 100 deg. There is a group of 
four faint apparent companion objects running North-South 
at a projected distance of ~ 40 kpc. 

2135-147 (PKS 2135-14, PHL 1657) 

The image of this source shown in Fig A9 shows a large 
companion galaxy at a projected distance of ~ 30 kpc ESE 
and an apparent close-in companion, or secondary nucleus, 
at a distance of ~ 10 kpc. A recent analysis of the spectrum 
of the close-in companion by Canalizo & Stockton (1997) 
has shown that this object is actually a foreground star. 
This object has also been imaged in the V-hand with HST 
by BKS. Their analysis also reveals the underlying host to 
be best matched by an elliptical galaxy model. This quasar 
was identified as an X-ray source in the ROSAT All-Sky 
Survey (Brinkmann et al. 1997). 

2355-082 (PKS 2355-082, PHL 6113) 
The new R-hand HST image of this source in Fig AlO reveals 
the presence of a apparently early-type host with a PA of 
w 180 deg. A group of five small companion objects can be 
seen to the NE. This object was identified as an X-ray source 
in the ROSAT All-Sky Survey (Brinkmann et al. 1997). 

A 1.3 The Radio- Quiet Quasars 
0052-^251 (PG 0052-^251) 

Although not very obvious in the grey-scale representation 
used in Fig All, the host galaxy of this quasar can be seen to 
have clear spiral structure in the raw 7?-band image. There 
are two spiral arms present to the East and West of the 
nucleus with the Eastern arm being more extended and ap- 
parently terminating in a companion object. The overlying 
contours for this image reveal it to have a highly-luminous 
nuclear component. 

This object has recently been imaged in the i?-band 
by McLeod & Rieke (1994) and in J-, H- and K-h&nd by 
Hutchings & Neff (1997). Hutchings & Neff performed an 
analysis of the numbers, magnitudes and colours of the com- 
panion objects of this source, covering an angular extent 
comparable to the full WF2 image. They identify some 22 
companion objects brighter than K = 19.6, the vast mar- 
jority of which have colours consistent with mature stellar- 
populations. Comparing their near-infrared images with pre- 
vious optical studies, Hutchings & Neff note that there is no 
evidence for tidal structure made up of old stars, but con- 
clude that the host is in a phase of secondary star formation 
possibly induced by interaction with one of the close group 
of companions. 0052-1-251 was also included in the V-band 
HST imaging programme of Bahcall, Kirhakos & Schneider, 
who list the host morphology as spiral and identify many of 
the knots seen in the Eastern arm with Hll regions. Interest- 
ingly Bahcall, Kirhakos & Schneider comment that the inner 
regions of the surface-brightness profile of 0052-1-251 are well 
matched by an r^^* law, consistent with the best-fit ellipti- 
cal host from the K-hand imaging of Taylor et al. (1996), 
and suggestive that there may be a strong bulge component 
to the host galaxy. This is of course precisely what we have 
found through the construction pf a combined disc-l-bulge 
model for the galaxy as detailed in Section 5. Despite the 
fact that the eye is drawn to the high surface-brightness spi- 



ral arms in the HST image, we find that in fact 70% of its 
total i?-band light is contributed by the spheroidal compo- 
nent. Finally, we note that this quasar has also been identi- 
fied as an X-ray source in the ROSAT All-Sky Survey (Yuan 
et al. 1998). 

0204+292 

In the new _R-band image of this quasar the host galaxy 
appears to be elliptical with a well-defined position angle 
of « 90 deg. There are two companion objects on the sub- 
image shown in Fig A12, ~ 70 kpc to the WSW and ~ 45 to 
the NW respectively. This quasar was identified as an X-ray 
source in the ROSAT All-Sky Survey (Yuan et al. 1998). 

1549-^203 (LB 906, IE 15498-^203) 

It is immediately obvious from the new 7i-band image of this 
object that the host galaxy is small, and relatively faint com- 
pared to the nuclear component. Using a different grey-scale 
to that used in Fig A13, there is a suggestion of spiral-like 
features to the NW and SE of the nucleus, with the SE arm 
terminating at the apparent companion object which can 
be seen ~ 20 kpc to the E of the nucleus. A large, luminous 
elliptical galaxy can just be seen on the edge of the frame 
to the SW, with numerous fainter companions visible inside 
a radius of ~ 60 kpc. The full WF2 image would appear 
to show numerous companion objects although, as pointed 
out by Taylor et al. 1996, the density of the environment 
of this quasar is uncertain due to the presence of a nearby 
foreground cluster at z ~ 0.14 (Stocke et al. 1983). The first 
electronic images of this quasar were presented by Hutch- 
ings & Neff (1992). They imaged the object in both the V- 
and /-bands, detecting what looked like a bar structure run- 
ning North-South through the nucleus, while noting that the 
surface-brightness profile of 1549-1-203 was not exponential. 
This quasar was identified as an X-ray source in the ROSAT 
AU-Sky Survey (Yuan et al 1998). 

2215-037 (EX 2215-037) 

The image of this quasar shown in Fig A14 reveals the host 
to be compact and apparently undisturbed. A small appar- 
ently unresolved companion is detected only ~ 10 kpc to 
the SW of the quasar nucleus. A large elliptical galaxy can 
be seen some 65 kpc to the North with a second appar- 
ently unresolved companion ~ 50 kpc to the East. This ob- 
ject has been previously imaged on the HST PC using the 
F702W (wide R) filter by Disney et al. (1995). Using a two- 
dimensional cross-correlation modelling technique, Disney et 
al. found the host galaxy to be excellently matched by an 
elliptical galaxy model. Their wide-i? image confirms the ex- 
istence of the unresolved companion at ~ 2" separation from 
the nucleus and suggests that the environment of 2215—037 
resembles a poor cluster. This suggestion is supported by 
the image of the full WF2 chips which shows > 10 compan- 
ion objects inside a radius of ~ 100 kpc. This quasar was 
identified as an X-ray source in the ROSAT All-Sky Survey 
(Yuan et al. 1998). 
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Figure Al. The radio galaxy 0230-027 
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Figure A2. The radio galaxy 0307+169 
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Figure A3. The radio galaxy 1215-033 
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Figure A4. The radio galaxy 1215+013 
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Figure A5. The radio galaxy 1330+022 
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Figure A6. The radio galaxy 1342-016 
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Figure A7. The radio-loud quasar 1020-103 
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Figure A8. The radio-loud quasar 1217+023 
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Figure A9. The radio-loud quasar 2135—147 
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Figure AlO. The radio-loud quasar 2355-082 
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Figure All. The radio-quiet quasar 0052+251 
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Figure A12. The radio-quiet quasar 0204+292 
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Figure A13. The radio-quiet quasar 1549+203 
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Figure A14. The radio-quiet quasar 2215—037 
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APPENDIX B: LUMINOSITY PROFILES 

In this appendix we provide the observed and best-fitting 
model luminosity profiles for 14 AGN for which the new HST 
data are reported in this paper. The profiles are followed 
out to a radius of 10", which is representative of the typical 
outer radius used in the image modelling. In each figure 
the azimuthally averaged data arc indicated by circles, the 
azimuthal average of the best-fitting 2-dimensional model is 
shown by the solid line, and the dotted line indicates the 
contribution made to the surface-brightness profile by the 
point-source component of the model (after convolution with 
the PSF). 
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Figure B2. 



The Radio Galaxy 1215-033 
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Figure B3. 
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The Radio Galaxy 1215+013 
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Figure B4. 
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The Radio Galaxy 1330+022 
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Figure B8. 



The Radio Galaxy 1342-016 
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Figure B9. 
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Figure Bll. 



The Radio Quiet Quasar 2215-037 
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